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Abstract 

We discuss recent developments in the hydro-dynamic description of strongly cou- 
pled conformal field theories using the AdS/CFT correspondence. In particular, we 
review aspects of the fluid-gravity correspondence which provides a map between a 
class of inhomogeneous, dynamical, black hole solutions in asymptotically AdS space- 
times and arbitrary fluid flows in the strongly interacting boundary field theory. We 
explain how the geometric duals to the fluid dynamics are constructed in a boundary 
derivative expansion and use the construction to extract the hydrodynamic transport 
coefficients. In addition we also describe the recent developments extending the cor- 
respondence to incorporate matter fields and to non-relativistic systems. Based on 
lectures given at the CERN Winter School on Supergravity, Strings and Gauge Theo- 
ries, Geneva, Switzerland (February 2009). 
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1 Introduction 

One of the important questions in modern theoretical physics involves understanding the 
dynamics of strongly coupled quantum field theories. Not only is this of theoretical interest, 
but there is a large class of real world physical systems where conventional perturbation 
theory is a poor description of the actual physics. A case in point which is partly relevant 
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to the current discussion is the fascinating state of matter discovered in heavy ion collisions, 
the Quark Gluon Plasma (QGP), which is known to behave as a nearly ideal fluid. 

An important tool in the theoretician's toolkit to address strong coupling dynamics is 
the AdS/CFT correspondence [1, 2, 3], which provides us with a holographic reformulation 
of field theory dynamics in terms of classical gravitational dynamics (in a higher dimensional 
spacetime) . In general the AdS / CFT correspondence relates two seemingly disparate systems 
- it provides a deep inter-connection between an interacting quantum field theory on the one 
hand and string theory in a curved background on the other. The field theories of interest are 
typically large N gauge theories (obeying large N factorization in the planar limit) and the 
holographic dual is in terms of string theory in an asymptotically AdS background. While at 
generic values in parameter space one is dealing with two intrinsically complicated theories, 
at corners of parameter/coupling space one or the other description simplifies. Clearly, one 
has a simple field theoretic description when the coupling parameter is taken to be small; 
perturbation theory becomes reliable. In this regime the holographic dual description is 
in terms of string propagating in a highly curved background. In the opposite limit when 
we dial the field theory coupling to be strong, we simplify the string theory into classical 
(super-)gravity. 

There is a way of thinking about the AdS/CFT correspondence that chimes well with 
intuition for large iV gauge theories which is worth bearing in mind. On general grounds one 
expects that the large N limit of a quantum gauge theory behaves effectively classically i.e., 
quantum fluctuations are suppressed by 1/N. Said differently, the full quantum dynamics 
in this planar limit should be encoded in terms of a classical gauge field configurations, 
the "Master field" which controls the dynamics in this regime. While we have no concrete 
candidate for this Master field in large N QCD, for a wide class of supersymmetric gauge 
theories which arise as the world- volume theories on D-branes, the AdS / CFT correspondence 
identifies a candidate Master field - this is just string theory (or classical gravity) in an 
asymptotically AdS spacetime. In the strong coupling regime of the field theory the dynamics 
of single trace operators is captured completely by classical Einstein gravity coupled to other 
fields. 

In the course of these lectures we are going to be interested in a specific limit of the 
correspondence - we wish in particular to simplify the dynamics of the field theory to that of 
an effective classical fluid dynamics. As we shall see there is a precise sense in which this can 
be done for any interacting quantum field theory, by focussing on near-equilibrium dynamics 
and restricting attention to long wavelength physics. Under the holographic map we are 
led to consider a particular class of gravitational solutions, which turn out to be dynamical 
black hole spacetimes. This limit of the AdS / CFT correspondence which provides a concrete 
relation between the physics of fluids and that of gravity is what we call the fluid-gravity 
correspondence [4]. In the course of these lectures we will derive this correspondence and 
see its utility in various contexts. Before delving into the detail however it is worthwhile to 
pause and take stock of the reasons for why this is an interesting endeavour. 
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Fluid dynamics: As a classical dynamical system fluid dynamics provides interesting the- 
oretical challenges. It is well known that for non-relativistic incompressible viscous fluids 
described by the Navier-Stokes equations the issue of finding globally regular solutions re- 
mains a open challenge, see for e.g., [5]. At the same time fluid dynamical evolution shows 
very intriguing physics such as turbulence whose detailed understanding is still lacking. Fur- 
thermore, the behaviour of energy cascade in turbulent flows and the corresponding inverse 
cascade in lower dimensions are intriguing phenomena that beg for a better explanation. 

A holographic mapping of the fluid dynamical system into classical gravitational dynamics 
could in principle help in unearthing some of these mysteries, at the very least by providing 
a new perspective on the problem. To be sure, much of the physics of turbulence and global 
regularity are of interest in the context of non-relativistic, incompressible Navier-Stokes 
equations, while a natural realization of hydrodynamics in the fluid-gravity correspondence 
leads to relativistic conformal fluids. This however is not a primary obstacle, for as we will 
discuss towards the end, generalizations of the fluid-gravity correspondence to relax these 
constraints already exist. An obvious fantasy would be to hope that one can formulate a 
holographic dual of turbulence, but this is beyond the scope of these lectures. 

Gravitational solutions: Over the past decade new remarkable stationary solutions to 
higher dimensional gravity have been discovered and have served to highlight the limita- 
tions of the folk-lore about gravitational dynamics built on our intuition for gravity in four 
dimensions (for a review see [6]). An excellent illustration of this are the black hole unique- 
ness theorems which fail to generalize straightforwardly to higher dimensions. The problem 
of stability of the solutions in higher dimensions is also reveals some differences from the 
lower dimensional analogs. One encounters for instance Gregory-Laflamme instabilities for 
black strings and black branes [7] and closely related instabilities for spinning black holes 
[8]. Understanding the stability domain of a given classical solution is important to get a 
clear picture of the classical phase space of higher dimensional solutions. 

Fluid dynamics provides an interesting window to understand the physics of black holes. 
The idea of applying hydrodynamic intuition to black holes dates back to the works on the 
membrane paradigm [9, 10], wherein one modeled the black hole horizon by a membrane 
equipped with fluid like properties. More recently, analog models for black hole stability 
problem have been proposed whereby the Rayleigh-Plateau instabilities of liquid droplets 
was associated with Gregory-Laflamme instabilities of black holes [11]. However, in these 
applications the fluid dynamics is merely an analogy, a mnemonic to understand the qualita- 
tive physical details in the complex gravitational setting by invoking a simpler fluid model. 

The fluid-gravity correspondence however provides a real duality between the hydrody- 
namic description and the gravitational dynamics. This in particular implies that one can 
draw a precise quantitative connection between the two and thus enables us to understand 
aspects of the phase structure of black holes solutions and their stability in terms of the fluid 
model. More pertinently for our current discussion this holographic duality also allows us to 
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systematically construct dynamical black hole solutions. As we will see in the course of the 
lectures every fluid flow in the boundary field theory will map to a black hole spacetime in 
the bulk geometric description with a regular event horizon. In fact, the fluid-gravity corre- 
spondence will enable us to algorithmically construct black hole geometries given solutions 
to the fluid equations of motion. 1 

Relevance to real world physics: Theoretical understanding of the state of matter 
produced in heavy-ion collisions at RHIC (and perhaps soon at LHC), the QGP, requires 
knowledge of dynamics in strongly coupled QCD. Current understanding is that subsequent 
to the collision of the ions, the resulting constituents of the system rapidly thermalize and 
comes into local thermal equilibrium and thenceforth evolve according to hydrodynamics 
until the local temperature falls back below the deconfinement temperature and the QGP 
hadronizes. The hydrodynamic regime is characterized by a set of transport coefficients; 
in particular, since much of the flow in the QGP is the shear driven elliptic flow, it is the 
shear viscosity of the plasma that is of most relevance. Ideally, one would like to be able 
to start from the microscopic description in terms of QCD and be able to compute these 
transport coefficients. However, with the QCD coupling constant remaining strong near the 
deconfinement temperature one needs to find a way to go beyond perturbation theory. The 
obvious choice, Lattice QCD is somewhat handicapped in this respect since it is ill-equipped 
to deal with real time physics and Lorentzian correlation functions. 2 

The AdS/CFT correspondence provides a theoretical framework to understand some of 
the qualitative features of the hydrodynamics seen in QGP, by providing us with an efficient 
way to access strongly coupled physics in a class of superconformal field theories. While 
these field theories are qualitatively different from QCD in their vacuum, one might argue 
that at finite temperature some of these differences are perhaps mitigated. An interesting 
observation based on lattice simulations of QCD free energy, is that there appears to be 
a range of temperatures say between 2T C and 5T C (recall that T c ~ 175 MeV for QCD) 
where the energy density as a function of temperature shows Stefan-Boltzmann scaling with 
a numerical pre-factor which is about 80% of the free field value. This could be taken as 
prima-facie evidence for an effective description as a strongly coupled CFT for it is very 
similar to the situation in M = 4 Super- Yang Mills (SYM) theory, whose strong coupling 
free energy is exactly 3/4 the free field value [18]. However, other observables such as pressure 
deviate from the value predicted by conformal invariance, thereby weakening the analogy. It 
is therefore worth keeping in mind that the theories one is discussing are not quite QCD. 

That said it is rather remarkable that the only class of strong coupling field theories 
for which we can compute hydrodynamic transport coefficients exhibit a remarkable quan- 

1 In fact there are several results in the literature exploring the phase structure [12, 13, 14] and stability 
[15, 16] of black holes using the dual fluid dynamics in the context of field theories compactified on spatial 
circles with supcrsymmctry breaking boundary conditions. These boundary conditions break conformal 
invariance and the field theories are confining in the infra-red. 

2 For recent developments on lattice computations see [17] and references therein. 
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titative agreement with those arising from numerical fits to RHIC data. For instance the 
universal behaviour of the shear viscosity in hydrodynamic description of field theories with 
gravitational holographic duals [19] has already attracted attention and has impacted exper- 
imental analysis of RHIC data, see for example [20, 21, 22, 23]. 3 In any event, independent 
of applications to heavy-ion collisions one can view the superconformal field theories as toy 
models; it is certainly quite remarkable that the holographic map allows us to explicitly 
determine the transport properties of a strongly coupled non-abelian plasma. 

Summary of the lectures: In these lectures we use the AdS/CFT correspondence to 
study the effective description of strongly coupled conformal field theories at long wave- 
lengths. On physical grounds it is reasonable that any interacting quantum field theory 
equilibrates locally at high enough energy densities, and so admits an effective description 
in terms of fluid dynamics. The variables of such a description are the local densities of 
all conserved charges together with the local fluid velocities. The equations of fluid dynam- 
ics are simply the equations of local conservation of the corresponding charge currents and 
energy-momentum tensor, supplemented by constitutive relations that express these currents 
as functions of fluid mechanical variables. As fluid dynamics is a long wavelength effective 
theory, these constitutive relations are usually specified in a derivative expansion. At any 
given order, thermodynamics plus symmetries determine the form of this expansion up to a 
finite number of undetermined coefficients. These coefficients may then be obtained either 
from measurements or from microscopic computations. 

The best understood examples of the AdS/CFT correspondence relate the strongly cou- 
pled dynamics of certain (super-) conformal field theories to the dynamics of gravitational 
systems in AdS spaces. In particular, we will demonstrate that Einstein's equations with a 
negative cosmological constant, supplemented with appropriate regularity restrictions and 
boundary conditions, reduce to the nonlinear equations of fluid dynamics in an appropriate 
regime of parameters. We provide a systematic framework to construct this universal non- 
linear fluid dynamics, order by order in a boundary derivative expansion i.e., as an effective 
theory. 

There is a rather rich history of studying hydrodynamics of non-abelian plasmas using 
holographic methods provided by the AdS/CFT correspondence. The early work of [25] was 
the first to relate the process of thermalization in the field theory with the study of black hole 
quasi-normal modes. Subsequently, the seminal work of Policastro, Son and Starinets [26] 
began the investigation of linearized fluid dynamics from linearized gravity in asymptotically 
AdS black hole backgrounds. The exploration of linearized hydrodynamics has been carried 
out in various different contexts over the years in [19, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 
37, 38, 39, 40, 41, 42, 43, 44, 45, 46] and we refer the reader to the excellent review [47] for 
other references on the subject. The fluid-gravity correspondence [4] itself was motivated in 

3 For a recent account of the nearly ideal fluids in encountered nature, such as the QCP and cold atoms 
at unitarity, see [24]. 
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part by these studies and the attempts to construct the holographic dual of the so called 
Bjorken flow [48, 49] 4 which is believed to be relevant to understanding the central region of 
heavy- ion collisions hydrodynamically [64]. At the same time the investigations of the fluid 
dynamical regime of stationary black holes in asymptotically AdS spacetimes [65] paved the 
way for a clear understanding of the hydrodynamic regime in the gravitational context. 

The fluid-gravity correspondence was originally discussed in [4] in the context of gravita- 
tional duals of four dimensional superconformal field theories whose holographic dual is given 
by AdSs xl 5 where X 5 a Sasaki-Einstein manifold which determines the CFT. A special 
case is the J\f = 4 SYM where X 5 = S 5 . In [66] the global aspects of the bulk geometry were 
discussed and a geometric construction of the fluid entropy current was given. Subsequently, 
this discussion was generalized to other dimensions in [67, 68, 69]. In addition [70] described 
how to include external forcing in the hydrodynamic description by placing the fluid on a 
curved manifold (and also explicitly included dilaton couplings). There is also a discussion 
of charged fluid dynamics [71, 72, 73, 74], which in the context of M = 4 SYM corresponds 
to looking at the grand-canonical ensemble with chemical potential for U(1)r charges, and 
in addition inclusion of magnetic and dyonic charges in AdS4 [75, 76]. While all these dis- 
cussions are for conformal fluids whose duals are asymptotically AdS spacetimes, one can 
extend the discussion to non-conformal fluids living on Dp-brane world-volumes [77, 78] as 
well as explorations involving higher derivative terms in the gravitational description [79]. 
Finally, the restriction of relativistic invariance can be relaxed to consider non-relativistic 
conformal fluids [80] as well as duals to incompressible Navier-Stokes flow [81]. For a review 
of some of the developments on the subject see also [82]. We will mainly review the basic 
features of the correspondence in these lectures and will discuss some of the generalizations 
towards the end in §6. 

The essential physical points arising from the fluid-gravity correspondence which we will 
describe in detail below can be summarized as follows: 

• The gravitational derivation of the relativistic Navier-Stokes equations and its higher- 
order generalizations confirms the basic intuition that fluid dynamics is indeed the 
correct long-wavelength effective description of strongly coupled field theory dynamics. 

• The geometries dual to fluid dynamics turn out to be black hole spacetimes with 
regular event horizons [66]. This indicates that the hydrodynamic regime is special 
and in particular can be interpreted to indicate that the fluid dynamical stress tensors 
lead to regular gravity solutions respecting cosmic censorship. 

• The explicit construction of the fluid dynamical stress tensor leads to a precise deter- 
mination of higher order transport coefficients for the dual field theory. The results we 

4 For further explorations of the spacetime geometry dual to Bjorken flow see [50, 51, 52, 53, 54, 55, 56, 
57, 58, 59] and [60] for a review. More recently, these class of geometries have been understood within the 
framework of the fluid-gravity correspondence in [61, 62, 63]. 
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present for the transport coefficients for four dimensional superconformal field theories 
with holographic duals were also derived in [83]. 

Outline of the lectures: We will begin with a review of relativistic fluid dynamics in §2 
and then discuss various aspects of conformally invariant fluids in §3, reviewing in particular 
the extremely useful Weyl covariant formalism developed in [84] in §3.2. In §4 we will review 
the basic scheme to construct gravitational solutions dual to fluid dynamics following [4]. 
We then turn to a discussion of the physical properties of our solutions in §5. Finally in §6 
we will discuss various generalizations of the fluid-gravity correspondence and conclude with 
a discussion in §7. 

2 Elements of fluid dynamics 

Fluid dynamics is the low energy effective description of any interacting quantum field theory, 
valid for fluctuations that are of sufficiently long wavelength. This description is intrinsi- 
cally statistical in nature, for it is the collective physics of a large number of microscopic 
constituents. Usually one thinks of the hydrodynamic description as the effective continuum 
model valid on macroscopic scales. For an excellent introduction to the basic ideas in fluid 
dynamics see [85] (see also [86] for relativistic fluids). 

To understand the statistical origins of the fluid dynamical description, consider a quan- 
tum system in a grand canonical ensemble, where we prescribe the temperature and chemical 
potentials for various conserved charges. In global thermal equilibrium we calculate observ- 
ables by computing correlation functions in the grand canonical density matrix. One can also 
perturb away from this global equilibrium scenario and allow the thermodynamic variables 
to fluctuate. For fluctuations whose wavelengths are large compared to the scale set by the 
local energy density or temperature, one describes the system in terms of fluid dynamics. 

One can pictorially think of the situation as follows; sufficiently long- wavelength fluctu- 
ations variations are slow on the scale of the local energy density/temperature. Then about 
any given point in the system we expect to encounter a domain where the local temperature 
is roughly constant - in this domain we can use the grand canonical ensemble to extract 
the physical characteristics of the field theory. Of course, different domains will described 
by different values for the intrinsic thermodynamic variables. Fluid dynamics describes how 
these different domains interact and exchange thermodynamic quantities. 

A more formal way to define the hydrodynamic regime is the following. In any interact- 
ing system there is an intrinsic length scale, the mean free path length £ m f p , which is the 
characteristic length scale of the interacting system. This is most familiar from the kinetic 
theory picture of gases, but applies equally well to fluids. In the kinetic theory context the 
£ m f p simply characterizes the length scale for the free motion of the constituents between 
successive interactions. To achieve the hydrodynamic limit we are simply requiring that 
we examine the system at length scales which are large compared to the £ m { p , so that the 
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microscopic inhomogeneities are sufficiently smeared out. 

The dynamical content of the hydrodynamic equation is simply conservation of energy 
and other conserved global charges. These can be succinctly summarized by giving the stress 
tensor T^, which is a symmetric two tensor, and the charge currents Jj , where I = {1, 2, • ■ ■ } 
indexes the set of conserved charges characterizing the system. The dynamical equations 
then are given as 



To specify the system further we just need to find expressions for the stress tensor and 
charge currents. Since the fluid dynamical description is only valid when the underlying 
microscopic QFT achieves local thermal equilibrium, we should be able to use the basic 
thermodynamic variables to characterize the system. Furthermore, we need to describe 
how the different domains of local thermally equilibrated fluid interact with each other. To 
understand how this can be achieved, let us focus on a fluid element which is characterized 
by the local values of the thermodynamic variables. As the fluid element can exchange its 
characteristics with neighbouring fluid elements, we should associate to this fluid element a 
velocity field u^, which describes the flow of thermodynamic quantities. It turns out that 
the thermodynamic variables together with the velocity field serve to characterize the fluid 
completely. 

Consider a QFT living on a <i-spacetime dimensional background 3d with (non-dynamical) 
metric g^ u . The coordinates on Ed will be denoted henceforth as x^. We can summarize the 
dynamical degrees of freedom characterizing the hydrodynamic description of this interacting 
QFT as: 

• Extrinsic quantities: Local energy density p and charge densities qi. 

• Fluid velocity w M (normalized g^ v u^u u = — 1). 

• Intrinsic quantities: Pressure P, temperature T, and chemical potentials /zj determined 
by equation of state. 

All that remains is to express the stress tensor and charge currents in terms of these variables. 
2.1 Ideal fluids 

Let us focus first on the description of an ideal fluid which has no dissipation. Then by 
passing to a local rest frame, where we choose the velocity field to be aligned in the direction 
of energy flow, we can identify the components of the stress tensor as the energy density 
(temporal component longitudinal to the flow) and pressure (spatial components transverse 
to the velocity field). Similarly, in this local rest frame the components of the charge current 
are the charge density itself (along the velocity field). Putting this together, we learn that 
for an ideal fluid 





. 





ideal 



pu^u u + P {gT + u»u u ) , 



(2.2) 
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Before proceeding to incorporate the physics of dissipation, let us pause to introduce 
some notation that will be useful in what follows. Since the ci-velocity w M is oriented along 
the temporal direction, we can use it to decompose the spacetime into spatial slices with 
induced metric 

P^ = g^ + u»u u . (2.3) 

We can view P^ v as a projector onto spatial directions; it is a symmetric positive definite 
tensor which satisfies the following identities: 

P*" u, = 0, P>» P pv = P% = P"> g pu , P/ = d - 1. (2.4) 

In terms of this projector we can express the ideal fluid stress tensor more simply as 

(Tn^ = P^u» + PP^. (2.5) 

While the decomposition of T^ u into the temporal and spatial components is not very in- 
teresting at the level of ideal fluids, it will play a more important role when we include the 
effects of dissipation. 

In addition to the stress tensor and charge currents, there is another quantity which we 
would like to keep track of, viz., the entropy current. Pictorially the entropy current keeps 
track of how the local entropy density varies in the fluid. For an ideal fluid given the entropy 
density s{x) the entropy current takes the simple form 

(^)idc al = ^- (2-6) 

Using the equations of motion (2.1) for an ideal fluid (2.2) and standard thermodynamic 
relations, it is easy to show that the entropy current is conserved 

V,(J£) ideal = 0. (2.7) 
i.e., the fluid flow involves no production of entropy. 

2.2 Dissipative fluids 

Ideal fluids with stress tensor given by (2.5) are an approximation; they don't include any 
physics of dissipation. This is clearly seen by the conservation the entropy current. In general 
we expect that the flow of the fluid results in the creation of entropy consistent with the sec- 
ond law of thermodynamics. More pertinently, dissipation is necessary for a fluid dynamical 
system to equilibrate when perturbed away from a given equilibrium configuration. Micro- 
scopically one can understand the dissipative effects as arising from the interaction terms 
in the underlying QFT. As a result the terms incorporating the effects of dissipation in the 
stress tensor and charge currents will depend on the coupling constants of the underlying 
QFT. 
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To model a hydrodynamical system incorporating the effects of dissipation we only need 
to add extra pieces to the stress tensor and charge currents. Let us denote the dissipative 
part of the stress tensor by IP" as the corresponding part of the charge currents by T M : 

(dissipative = ?/«" + T". (2.8) 

To complete specification of the hydrodynamic system we should determine IT M " and T M in 
terms of the dynamical variables the velocity field, u^, and the thermodynamic variables p, 
P, qi, etc.. 

The traditional way to determine the dependence of the dissipative components of the 
stress tensor and charge currents is somewhat phenomenological. One employs the second 
law of thermodynamics (positive divergence of the entropy current), cf., [85, 86], to determine 
the set of allowed terms in the most general form of the constitutive equations. While this 
makes some sense for incorporating the leading non-linearities in the fluid description, at 
higher orders the procedure starts to be fraught with ambiguities as exemplified by the Israel- 
Stewart formalism [87, 88]. We will take a different viewpoint on this and in fact outline a 
procedure to unambiguously construct a non-linear stress tensor for hydrodynamics. 

The procedure we have in mind is inspired by the manner one constructs effective field 
theories. When one writes down an effective field theory for a QFT, at any given order 
one takes into account all possible terms that can appear in the effective lagrangian con- 
sistent with the underlying symmetry. These irrelevant operators are suppressed by power 
counting in momenta and appear with arbitrary coefficients in the effective lagrangian. We 
have argued that hydrodynamics is best thought of as an effective classical description of 
any interacting field theory valid when the system achieves local thermal equilibrium. Hence 
following the logic of effective field theories we should allow in the hydrodynamic approxima- 
tion all possible operators consistent with the symmetries. These irrelevant operators in the 
present case are nothing but derivatives of the velocity field and thermodynamic variables. 
Thus one can cleanly express the fluid dynamical stress tensor in a gradient expansion. 5 

Before proceeding to enumerate the various operators that can contribute to the stress 
tensor, let us constrain the dissipative components. First of all we have a choice of frame 
which is related to how we choose the fluid velocity. Secondly when we enumerate gradient 
terms at any given order we are free to use the lower order equations of motion to reduce 
the set of gradient terms. Let us discuss these two issues in turn and see how they serve to 
constrain the dissipative terms. 

Defining the velocity field: Since we are discussing relativistic fluids, the fluid velocity 
field needs to be defined with some care. The issue is simply that in relativistic dynamics one 

5 The power counting scheme we have in mind for fluid dynamical operators is the following: each spacetime 
derivative will count as dimension 1 for book keeping purposes. It is important to bear in mind that the 
actual scaling dimension of the operator (as discussed in §3) which follows from the transformation of the 
stress tensor under scale transformations is different. 
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cannot distinguish between mass and energy in a clean fashion. Usually in non-relativistic 
fluids one would talk about heat flow between different fluid elements, but for a relativistic 
fluid heat flux necessarily leads to mass or momentum transfer and hence to energy flow. 
By picking the velocity field appropriately we can fix this ambiguity; the precise fashion in 
which we choose to do so is simply a matter of convention. 

For ideal fluids we picked the velocity field such that in the local rest frame of a fluid 
element, the stress tensor components longitudinal to the velocity gave the local energy 
density in the fluid. One can in fact demand the same to be true when we incorporate the 
effects of dissipation. This leads to the so called Landau frame which is defined by demanding 
that the dissipative contributions be orthogonal to the velocity field i.e,. 



Formally, one can define the Landau frame as follows: the stress tensor which is a symmetric 
two-tensor on the background manifold Bd has a single timelike eigenvector. In the Landau 
frame we define the velocity field m m to be given by this eigenvector (unit normalized), so 
that the definition of the velocity field is tied to the energy-momentum transport in the 
system. For charged fluids one has the possibility of working in the so called Eckart frame 
where the fluid velocity is tuned to charge transport - here we define the velocity field to 
be the determined by the charge current. In what follows we will work exclusively in the 
Landau frame and hence apply (2.9) to constrain the dissipative terms of the stress tensor 
and charge currents. 

Enumerating independent operators: The issue of enumerating the independent op- 
erators at any given order is completely analogous to the discussion in effective field theory. 
Suppose we are interested in generalizing the ideal fluid (2.2) to first order in the gradient 
expansion. For the stress tensor we would write down all possible symmetric two tensors 
built out of the gradients of the velocity field and thermodynamic variables. However, the 
ideal fluid equations of motion (2.1) themselves are first order in derivatives, and therefore 
relate the derivatives of the thermodynamic potentials to the gradients of the velocity field. 
Employing these relations we can simplify the expression for the first order stress tensor 
to be given in terms of derivatives of the velocity field alone. This process can clearly be 
iterated to higher orders. 

Decomposing velocity gradients: We are interested in analyzing various combinations 
of derivatives of the velocity field. First let us consider the decomposition of the velocity 
gradient V v u^ into a part along the velocity, given by the acceleration a M , and a transverse 
part. The latter can in turn be decomposed into a trace, the divergence 8, and the remaining 
traceless part with symmetric and antisymmetric components, respectively given by the shear 
a^ u and vorticity u^. The decomposition can then be written as follows: 




(2.9) 



Vu" = 



1 



(2.10) 



d- 1 
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where the divergence, acceleration, shear, and vorticity, are defined as: 

e = vx = pliV 

r \ t \ 1 R 1 (2.11) 

= V^u v) + u {fl a v) 9 = P^ P vP V (a u 0) 9 P" u K J 

a — 1 a — 1 

a"" = VK" ] + u [ " a u] = P» a P vP V [aU(3] . 

For future reference we note that we will also have occasion to use a the following notation to 
indicate symmetric traceless projections transverse to the velocity field. For any two tensor 
T^ v we define: 

T (H = P m pvP T{ap) _ — P ^ pa/3 Tafj _ ^ ^ 

Note that we can write the projectors a bit more compactly: P^ a P v @ V( a = P^'SJ p u u > 
and P^ a P uP V [a u m = P*V ' p u u l It is easy to verify all the previously asserted properties, 
in addition to a? = and P^ v = a u : 

We are now in possession of sufficient amount of data to write down the dissipative part 
of the stress tensor to leading order in the derivative expansion. First of all let us notice 
that the zeroth order equations of motion i.e., those arising from the ideal fluid description 
relate the gradients of the energy density and pressure to those of the u^. The quickest 
way to derive the required relation is to consider projections of the conservation equation 
V M (P M ^) idcal = 0, along the velocity field and transverse to it, i.e., 

^V M (Tn idcal = =► (p + P)V/ + ^ = 
4V,(r) ideal = CV/+(/) + F)P ra /V/ = 0. (2.14) 

respectively. To characterize the stress-tensor at leading order in the gradient expansion 
our task is reduced to writing down symmetric two tensors that can be built solely from 
velocity gradients and we should furthermore account for the Landau frame condition. These 
conditions in fact isolate just two terms which can appear in the expression for IP 1 ': 

= -2 rj - C 9 , (2.15) 

where we have introduced two new parameters the shear viscosity, r\, and the bulk viscosity, 

c 

Likewise for the charge current T M we will obtain contributions which are first order in the 
derivatives of the thermodynamic variables p and qi and also the velocity field (where now 



6 Note that we use standard symmetrization and anti-symmetrization conventions. For any tensor F a f, 
we define the symmetric part Fu^ = \ (F a j, + Fj, a ) and the anti-symmetric part Fi,m = \ {Fab — Fba) 
respectively. We also use & to indicate the velocity projected covariant derivative: = it M V M . 
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the acceleration contributes). We will choose this time to eliminate the acceleration term 
using zeroth order equation of motion (2.14) and express derivatives of energy density and 
charge density. Usually these are the only contributions which are considered for the charge 
current. However, there is another potential contribution from a pseudo- vector [71, 72], 

f = ^«"VV, (2.16) 

provided we are willing to mix contributions with different parity structure at the first order. 
In fact, for fluid dynamics derived from the AdS/CFT correspondence for the M = 4 Super 
Yang-Mills fluid, there is indeed a contribution of this sort to the charge current (in a specific 
ensemble). 7 Putting things together we find that we have a first order contribution to the 
charge current: 

Tfo = -k 13 V u q,j - 7/ V»P ~ 0/ P , (2.17) 

where we again have new coefficients, xjj is the matrix of charge diffusion coefficients, 7/ 
indicates the contribution of the energy density to the charge current, and 0/ which are the 
pseudo-vector transport coefficients. For later purposes, it will be convenient to re-express 
the charge current in terms of the intensive parameters \xi and T (chemical potential and 
temperature), so that (2.18) can be recast equivalently as 8 

Tf 1}J = -yen P» v V„ (^Q - U T - 7/ P^ V„T . (2.18) 

Assembling all the pieces together we claim that a generic charged fluid flow will satisfy 
the dynamical equations (2.1) with the stress tensor and charge currents at leading order in 
gradient expansion (dropping therefore higher derivative operators): 

= pu^vT + P (g^ + u>*u v )-2ria» v -(9P» V , 
J* = qm^-xuP^Vv (|Q-0/^-7/^V,T. (2.19) 

Following conventional terminology we will refer to the fluid dynamical system specified by 
(2.19) as a viscous fluid or relativistic Navier-Stokes equations. 9 

Apart from the thermodynamic potentials, we have a set of transport coefficients 77, £, 
x/j, 7/ and U 1 which need to be determined to completely specify the relativistic viscous 
fluid. If the QFT whose hydrodynamic description we seek is weakly coupled then we can 
in principle determine these coefficients in perturbation theory. This is per se not a trivial 
exercise, for a discussion of the computation of transport coefficients in perturbative gauge 
theory see [89, 90]. We will however be interested in deriving these coefficients for a strongly 
coupled quantum field theory using a dual holographic description. 

7 This contribution can be traced to Chcrn-Simons couplings of the bulk gauge field in the gravitational 
description, c.f.. [7f, 72]. 

8 One can relate the transport coefficients {xjj,jj} to the original coefficients {kjj^jj} in terms of 
various susceptibilities. However, getting explicit expressions for these in general will require knowledge of 
the precise thermodynamic functions. We thank R. Loganayagam for useful discussion in this regard. 

9 We will review the conventional non-relativistic equations later, see §6.2. 
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Finally, let us turn to the other quantity of interest, viz., the entropy current. Including 
the effects of dissipation we learn that entropy is no longer conserved. Under fluid evo- 
lution we should have entropy production since the system is generically evolving from a 
non-equilibrium state to an equilibrium state. Entropy increase consistent with the second 
law can be rephrased as the statement that the entropy current should have non-negative 
divergence. For a relativistic viscous fluid we will therefore have an entropy current Jg with 

J£ > . (2.20) 

It is possible to show that for uncharged fluids involving pure energy momentum transport 
one has a entropy current statisfying: 

V,r s = ^a aP a^ , (2.21) 

which is non-negative as long as the shear viscosity is positive. We will return to a detailed 
discussion of the entropy current for conformal fluids in §5.4. 

2.3 Causality issues in relativistic viscous fluids 

A natural question is whether the relativistic viscous fluids can be extended to higher or- 
ders. 10 To determine this we should just follow the conventional philosophy of effective field 
theory and try to write down the full set of operators at the desired order in gradient ex- 
pansion. These higher order terms will be less and less important (they are irrelevant in the 
low energy effective theory) as we go to longer wavelengths and are thus somewhat inconse- 
quential for the hydrodynamic evolution. The leading order terms which we incorporated in 
the viscous fluid (2.19) are singled out because of their importance in realizing a channel for 
the fluid to relax back to equilibrium - they are the leading dissipative terms in the theory. 

However, in conventional relativistic fluid literature one usually encounters a different 
rationale having to do with causality issues and the initial value problem. Consider the 
viscous fluid system (2.19) with the dynamical equations (2.1). This system of equations is 
first order in time derivatives, 11 and as a result the system of partial differential equations 
equations describing a relativistic viscous fluid is parabolic. 

In order to have a relativistic system with a well defined initial value problem the equa- 
tions of motion are required to be hyperbolic. Only for hyperbolic partial differential equa- 
tions is it possible to evolve initial data specified on a Cauchy surface. Usually lack of 
hyperbolicity leads to formation of singularities and manifests itself by acausal behaviour. A 
useful analogy to keep in mind, which is in fact central to the current discussion, is the heat 
equation, which being first order in time derivatives is parabolic. Solutions to the heat equa- 
tion show diffusive behaviour which leads to acausal propagation of signals (when viewed 
with respect to the light-cone of the underlying Lorentzian manifold Bd). 

10 We will refer the generalized system of fluid dynamical equations as non-linear viscous hydrodynamics, 
notwithstanding the fact that the relativistic Navicr-Stokes equations themselves are non-linear. 

11 The simplest way to see this is to work in the local inertial frame and realize that the gradient terms in 
(2.19) only involve spatial derivatives owing to the projector P^ v ' . 
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This issue has led several people to come up with so called "causal relativistic hydro- 
dynamics", the most prominent among which is the phenomenologically motivated Israel- 
Muller- Stewart formalism [87, 88]. 12 For sake of completeness we will briefly review the 
Israel-Stewart formalism and in particular will argue that the issue of causality violation in 
the hydrodynamic description is a red-herring. 

The Israel- Stewart idea is to add higher order terms to the relativistic viscous fluid system 
(2.19). The terms that are added are such that the resulting entropy current satisfies the 
second law (2.20). However, it is worth noting that their phenomenological formalism doesn't 
fully account for all the irrelevant operators one can add at the corresponding order [83] . We 
will illustrate the formalism with a simple toy model which captures the basic idea behind 
their approach. 

Consider the derivation of the diffusion equation from the continuity equation, which is 
a conservation equation 13 

^ + V-j = 0, (2.22) 

and the phenomenological Fick's law 

J = -DVp, (2.23) 
where D is the diffusion constant. Eliminating the current j we find the diffusion equation 

^- J DV 2 p = 0. (2.24) 

The equation (2.24) of course has a diffusive spectrum uo ~ —ik-k and is a simple example 
of a parabolic PDE which doesn't have a well defined initial value problem. The problem is 
that any perturbation to the system instantaneously spreads out to arbitrary spatial scales - 
essentially the fluctuations don't relax quickly enough. In the diffusion equation (2.24), due 
to absence of quadratic term in time derivatives, the velocity of propagation of disturbances 
is infinite. Equivalently, in the absence of density fluctuations there is nothing forcing the 
current to relax to its stationary equilibrium value. Ideally, the current should damp out 
exponentially in time. This can be achieved by phenomenologically modifying Fick's law 
(2.23). Consider introducing a relaxation time scale T n so that (2.23) is modified to 

3 = -DVp- r w d t j , (2.25) 

One can think of the deformation to Fick's law described above as an example of the leading 
irrelevant operator that we include in the system. Now eliminating the current between 
(2.22) and (2.25) we find a hyperbolic equation, for 

^-DV 2 p-r 7T V-d t j=0 ^ ^-DV 2 p + r n d 2 p = . (2.26) 



12 For reviews on this subject we refer the reader to [86, 91]. In the context of heavy-ion collisions an 
excellent discussion complemtenting our presentation here can be found in [92]. See also [58] for related 
discussion in the holographic context. 

13 For the sake of simplicity we will momentarily switch to a non-relativistic system; for the present 
discussion V is a spatial gradient. 
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Note that we have used the fact that we are interested in the lowest order terms in the 
gradient expansion and are therefore free to use the lower order equations of motion (2.24) 
to simplify the dynamical equations. (2.26) has a finite velocity for the propagation of the 
density fluctuations, v pvop = ^JUJt^. 

The toy model illustrates the basic idea involved in the Israel- Stewart construction of 
"causal relativistic hydrodynamics" - in their construction one adds a class of second deriva- 
tive terms to the viscous relativistic hydrodynamics (2.19). It transpires that their construc- 
tion actually ignores some of the two derivative operators allowed by the symmetries, see 
the discussion in [83] in the context of conformal fluids. 

However, as clearly described in [93] this whole discussion is rather misleading. The 
essential point is that the effects of the higher order terms which are added to render the 
theory causal are of decreasing relevance deep into the hydrodynamic regime. From the 
effective field theory viewpoint emphasized here this is completely obvious - the irrelevant 
operators have diminishing role to play as we move into the deep infra-red. Modes which 
manifest acausal behaviour are outside the long-wavelength regime and are not part of the 
hydrodynamic description. Curing this causality problem completely requires us to actually 
go back to the microscopic description of the theory. It is nevertheless interesting to ask 
what the systematic procedure to extend the relativistic viscous hydrodynamics is and we 
will address this question to second order in derivatives in the rest of these lectures. 

3 Conformal fluids 

So far our discussion of fluid dynamics has been quite general and is valid for the hydrody- 
namic limit of any interacting Lorentz invariant quantum field theory. We now turn to a 
discussion of field theories which are conformally invariant. After briefly reviewing the con- 
formal transformation properties of the stress tensor, we turn to the construction of operators 
in the gradient expansion that transform homogeneously under conformal transformations. 
This will allow us to classify the set of operators that can appear in conformal hydrodynam- 
ics beyond leading order; for pure energy momentum transport this was first undertaken in 
[83]. We will then undertake a brief overview of an extremely useful framework for discussing 
conformal fluids, the so called Weyl covariant formalism [84]. This will be useful when we 
construct gravitational duals to conformal fluids in §4. 

3.1 Weyl transformation of the stress tensor 

Consider a relativistic fluid on a background manifold Bd with metric g^ v . We wish to 
consider a local Weyl rescaling of the metric and understand the transformation of the stress 
tensor. Given a Weyl transformation of the boundary metric 

9»u = e 2 *g, v =► g*' = e- 3 *V a ' 1 (3.1) 
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it is clear that the velocity field appearing in the stress tenor transforms as 

v p = e -*u li , (3.2) 

which is a direct consequence of the normalization of the velocity field, u^u^ = — 1. It 
follows from (3.1) and (3.2) that the spatial projector transforms homogeneously as well: 

ptiv _ giiv _|_ u V u v — g-20 pV-v 

We are interested in the behaviour of the stress tensor under conformal transformations. 
In general a tensor Q with components Q^'-'-u^ is said to be conformally invariant if it 
transforms homogeneously under Weyl rescalings of the metric, i.e., Q = e~ w ^ Q under (3.1). 
The real number w is the conformal weight of the tensor. It is important to remember that 
the weight of a tensor operator under Weyl transformations depends on the index positions. 14 
We also should require that the dynamical equations satisfied by Q remain invariant under 
conformal transformations. 

It is easy to work out the constraints on the stress tensor using the dynamical equation 
at hand (2.1). We find that first of all the conformal invariance requires that the stress 
tensor be traceless = and that it transforms homogeneously under Weyl rescalings of 
the metric with weight d + 2 

pfiv _ e -(d+2)4> pflU (3 3) 

These statements follow from the conservation equation (2.1) for the stress tensor, written 
in the Weyl transformed frame. We refer the reader to Appendix D of [94] for a derivation of 
these results. The tracelessness condition in fact supplies the equation of state for conformal 
fluids. Going back to the ideal fluid (2.2), we find that 

T; = =► P=J- [P . (3.4) 

This relation between pressure and energy density fixes the speed of sound in conformal 
fluids to be a simple function of the spacetime dimension: c s = /g— p Likewise it is not 
hard to show the charge current transforms homogeneously with weight d under conformal 
transformations 

J* = e -^>. (3.5) 

Another piece of data we need is the information of the scaling dimensions of the ther- 
modynamic variables. It is for instance easy to see that the temperature scales under the 
conformal transformation with weight 1. Similarly it follows from the thermodynamic Gibbs- 
Duhem relation, P + p = sT + fii qj, that the chemical potentials fii also have weight 1. 
Hence using (3.3) and (3.1) it follows that the energy density is given by a simple scaling 
law (a natural extension of the Stefan-Boltzmann scaling to d spacetime dimensions) 

T = e~ 4> f p~T d (3.6) 



14 One could equivalently talk about the invariant conformal dimension which involves a shift by a simple 
linear combination of the number of upper and lower indices: w- mv = w — n\ owet + 7i U p per . This follows from 
accounting for the weight of the metric. 
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These results together imply that an ideal conformal fluid has a stress tensor 

{Tn^ = *T d {g^ + du»u v ) , (3.7) 

where a is a dimensionless normalization constant which depends on the underlying micro- 
scopic CFT. 

To construct the fluid stress tensor at higher orders we simply need to enumerate the 
set of operators in our gradient expansion which transform homogeneously under (3.1). Let 
us consider the situation at first order in derivatives explicitly. Using the fact that the 
Christoffel symbols transform as [94] 

r v = % + K + s; d x <f> - g x „ g»° , 
we can show that the covariant derivative of transforms inhomogeneously: 

VX = <9X + F ma « A = ^ [V„ v? + S; u° d a <p - g, x u x g v ° d a <p\ . (3.8) 

This equation can be used to derive the transformation of various quantities of interest in 
fluid dynamics, such as the acceleration a M , shear a^, etc.. 

9 = VX = e"* (VX + (d - 1) u a d a <^ = e" (d+ (d - 1) ^) , 

a" = $u v = u"VaU v = e~ 2 ^ (V + P va d a <f)) , 

V ' (3.9) 

a* ' = P x ^V x u u) - - P^ V x u x = e" 30 a"", 

a — 1 

where in the last equation we have accounted for the fact that all epsilon symbols should be 
treated as tensor densities in curved space. The objects with correct tensor transformation 
properties scale as metric determinants i.e., oc \/—g, and e al3lS oc 77= j from which 

it is easy to infer their scaling behaviour under conformal transformations; in particular, 
e a ^s = e aPlS and e a ^ s = e a ^ s . 

Armed with this information we can examine the terms appearing in the first order 
conformal fluid. From (3.9) we learn that the expansion 9, transforms inhomogeneously 
implying that the coefficient of the bulk viscosity should vanish for a conformal fluid ( = 
0. 15 Likewise, for the charge current we can show that the contribution from the chemical 
potential and temperature should be in the combination fii/T; the term involving gradient of 
the temperature P yw V U T transforms inhomogeneously under Weyl transformations requiring 
77 = 0. Therefore a conformal viscous fluid has a stress tensor which to first order in the 
gradient expansion takes the form 

T i»> = a T d (g^ + du^u u ) -2r]a^ , 
J? = qm^-xuP^Vvf^f) -VxP , (3.10) 
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Equivalcntly, this follows from the tracelessness of the stress tensor. 
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where we have used the generalized Stefan-Boltazmann expression for the pressure. We also 
can deduce the conformal properties of the transport coefficients; in particular, 

V = T d ^f)^j/T). (3.11) 

3.2 Weyl covariant formulation of conformal fluid dynamics 

It is straightforward to carry out this exercise to higher orders in the gradient expansion by 
again analyzing the transformation properties of the various operators under Weyl rescalings. 
This was carried out for uncharged fluids in [83] to second order and a detailed check of this 
for the specific second order stress tensor derived holographically is given in [4]. However, 
it is useful to work a little more abstractly and exploit a little more of the structure of the 
conformal symmetry to present the answer compactly in a Weyl covariant form. We therefore 
will pause to review some details of this presentation following the beautiful work of [84]. 

The basic point is as follows: when we are discussing a conformal field theory and the 
associated hydrodynamic description on a background manifold Bd we are not interested in 
the metric data on Bd, but rather on the conformal structure on this background manifold. 
We can exploit this fact to simplify the structure of the derivative expansion by directly 
working with the conformal class of metrics on Bd- For brevity, let us denote the background 
with this conformal class of metrics as (Bd,C). 

On (Bd, C) we will define a new derivative operator that keeps track of the Weyl trans- 
formation properties better. A key fact that we will exploit is that fluid dynamics comes 
equipped with a distinguished vector field, the velocity, which as emphasized earlier is just 
the timelike unit normalized eigenvector of the stress tensor. We start by defining a torsion- 
less connection called the Weyl connection, whose associated covariant derivative captures 
the fact that the metric transforms homogeneously under conformal transformations (with 
weight —2). In particular, the Weyl connection V Weyl requires that for every metric in the 
conformal class C there exists a connection one-form A^ such that 

V^ l g, v = 2A a g, u . (3.12) 

Given this derivative structure, we can go ahead and define a Weyl covariant derivative X> M 
as whose action on tensors transforming homogeneously with weight w under i.e., Q£." = 
e ~ w <P Qv - j s given by 

ZW: = v A Q£:: • wA x QZ::: 

+ ( 9 x a A^-5^A a -5^A x ) Q£"; H 

- {gx» A a -51 A v - 5: A x ) Q£: - ■ • • . (3.13) 

The nice thing about this definition is that T>\Q%\\\ = e~ w ^DxQv-.i i- e -; the Weyl covariant 
derivative of a conformally invariant tensor transforms homogeneously with the same weight 
as the tensor itself. 
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One can view the Weyl covariant derivative as being determined in terms of the Weyl 
connection via T)^ = V^ cyl + w A^. From this description it follows immediately that the 
Weyl connection is metric compatible T> a g^ u = which follows from (3.12) and the fact 
that w = —2 for g^ u . In addition, in fluid dynamics we will require that the Weyl covariant 
derivative of the fluid velocity be transverse and traceless, i.e., 

u a V a u^ = 0, V a u a = 0. (3.14) 

These conditions enable to uniquely determine the connection one-form A^ to be the the 
distinguished vector field 

= u x V\u^ - j- Up V x u\ = a M - j—[ 9u » ' ( 3 - 15 ) 

built from the velocity field of fluid dynamics. 

One can rewrite the various quantities appearing in the gradient expansion of the stress 
tensor in this Weyl covariant notation. For instance at first order in derivatives we have the 
objects constructed from the velocity field 

= V^u v) , = -VK" ] , (3.16) 

both of which can be seen to have weight w = 3. The dynamical equations of fluid dynamics 
(2.1) can be recast into the Weyl covariant form. For instance, stress tensor conservation 
can be easily seen to be £> M T^ u = 0; for 

= V M + wA, + ( 9im A" -S£A a - 6£ A,) T au + {g m A v -8^A a - 5 v a A,) 
= V M T^ + (w - d - 2) A^T^ - A" Tjl 

= V M T^ , (3.17) 

where we used the conformal weight w = d + 2 of the stress tensor and the tracelessness 
condition to remove the inhomogeneous terms. For CFTs on curved manifolds one has the 
possibility of encountering a conformal anomaly in even spacetime dimensions; incorporating 
the trace anomaly W we can write the fluid dynamical equation as 

= + A» (T^ -W) = . (3.18) 

3.3 Non-linear conformal fluids 

We will now discuss the general conformal fluids to second order in derivative expansion using 
the Weyl covariant formalism, with an aim to cast fluid mechanics in manifestly conformal 
language. In order to achieve this we need to write down the set of two derivative operators 
that transform homogeneously under conformal transformations. These operators involve 
either two derivatives acting on the dynamical degrees of freedom or terms involving squares 
of first derivatives. We now proceed to enumerate these operators: 
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For operators built out solely from the fluid velocity, one has the following terms which 
explicitly have two derivatives 



In addition we can have terms which have the combinations: 



(3.19) 



o» a o ua = e-^a^a™ , u^u™ = e^u^ST , a" a u va = e^a^u™ , (3.20) 

which involve the squares of the first derivative operators. 

In addition to the operators built out of the velocity field we should also consider terms 
which involve the temperature T and various chemical potentials Recalling that both of 
these have weight 1 under conformal transformations, i.e., T = e~^T and fij = jlj, we 
can write down the relevant operators involving no more than two derivatives: 

^ (f ) = ^ (|) > ^ T = e ~^ f 

^ ^ V ~ v (3.21) 

V x V a (^>j = V X V„ (^Lj V x V a T = e^V x V a T . 

To complete the classification of the various tensors that can be constructed at the second 
derivative level, we need to study the curvature tensors that appear via the commutators of 
two covariant derivatives. These appear because the fluid couples to the background metric 
g^u on the manifold Bd and the curvature terms arise at second order in derivatives. These 
terms will be of particular importance in the holographic context; for example in deriving 
the fluid dynamical behaviour of M = 4 SYM on S 3 x R we encounter curvature terms. We 
review the definition of the curvature tensors in the Weyl covariant formalism in Appendix 
A. For the present we just record the symmetric traceless tensor involving two derivatives 
of the background metric on Bd, 

C muP u a u? = C muP u a u? (3.22) 
where C a p^ v is the Weyl tensor see (A. 8). 



3.4 The non-linear conformal stress tensor 

We are now in a position to discuss the non-linear hydrodynamics for conformal fluids. For 
the sake of brevity we will confine attention to fluids which have no conserved charges. This 
is equivalent to working in the canonical ensemble and focussing just on energy momentum 
transport. The results we derive here will be universal for a wide class of conformal fluids. 

Let us classify all the operators which are Weyl invariant at various orders in the derivative 
expansion. The results of the previous section §3.3 can be summarized as follows: at first two 
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orders in derivatives the set of symmetric traceless tensors which transform homogeneously 
under Weyl rescalings are given to be 16 

First order : o^ v 

Second order : = 2u a Vad^ , Xif = u a vP , 

Xf = 4 a a{ " a v \ , Xf = 2 a a< ^ w*^ , Xf = uj v \ 

where we have introduced a notation for the second derivative operators which will be useful 
to write compact expressions for the stress tensor below. Armed with this data we can 
immediately write down the general contribution to the stress tensor as: 

= r^Xf + ^Xf + A x Xf + A 2 Xf + A 3 Xf . (3.23) 

There are therefore have six transport coefficients rj, r n , k, Aj for i = {1,2,3}, which char- 
acterize the flow of a non-linear viscous fluid. 

For a fluid with holographic dual using the explicit result of the gravitation solution (5.1) 
(see §5.1) and the holographic stress tensor [95, 96] we find explicit values for the transport 
coefficients. In particular, for the M = 4 SYM fluid one has 17 [4, 83] 

V = ^~ {nT) => - = — , 
Q7T s Air 

r - = ^^' K = Vf (3 - 24) 

Ax = -5-, A 2 = ^, A 3 = 0. 

where in the first line we have used the standard entropy density for thermal M = 4 SYM 
at strong coupling s = ^ N 2 T 3 to exhibit the famous ratio of shear viscosity to entropy 
density [19]. 

The general analysis in AdS^+i of carried out in [68, 69] in fact allows us to write down 
the transport coefficients described above in arbitrary dimensions in nice closed form. 18 The 
transport coefficients for conformal fluids in d- dimensional boundary 3d are: 

V 16nG { " +1) \ d J s 4vr' 



d 

Tn ~ 4^T 



1 x T (1 

1 + — Harmonic — 
d \ d 



d rj 
K ~ 2tt (d - 2) T 



(3.25) 



A 1 = ^|, A 2 = -L Harmonic 1, A 3 = 



16 Various papers in the literature seem to use slightly different conventions for the normalization of the 
operators. We will for convenience present the results in the normalizations used initially in [83] . This is the 
source of the factors of 2 appearing in the definition of the tensors Ti, see [4] for a discussion. 

17 The result for generic JV = 1 supcrconformal field theories which are dual to gravity on AdSs x J 5 are 
given by simply replacing by the corresponding central charge of the SCFT. 

18 For d = 3 the general results were initially derived in [67]. See also [59] for determination of some of 
these coefficients in d = 3 and d = 6. 
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where Harmonic (x) is the harmonic number function. 19 
4 Non-linear fluid dynamics from gravity 

Our discussion of fluid dynamics has thus far been very general and in fact is valid for 
any interacting field theory As we have seen one can distill the information in this low 
energy effective description into the specification of a finite number of transport coefficients. 
The transport coefficients can in principle be determined from the microscopic field theory. 
However, if we are interested in systems which are intrinsically strongly coupled then we 
need to find a technique to extract the non-perturbative values of the transport coefficients. 
In the rest of these lectures we will focus on such strongly coupled scenarios. As is well 
known, for a class of such theories one can exploit the AdS/CFT correspondence to extract 
the hydrodynamical parameters and this will be our prime focus. 

Let us consider a d dimensional field theory on a background manifold Ed which is holo- 
graphically dual to string theory on an asympotically AdS^+i spacetime (which is perhaps 
warped with some internal compact manifold to be a critical string background). The pro- 
totype example which we should keep in mind is the classic duality between Type IIB string 
theory on AdSs xS 5 and the dynamics of M = 4 SYM. Following standard parlance we will 
refer to the field theory as living on the boundary of the asymptotically AdS spacetime. 20 
In the example of M = 4 SYM we have two dimensionless parameters in the field theory, 
the 't Hooft coupling A and the rank of the gauge group N (we will consider unitary gauge 
group SU(N)). In the limit of large N and large A the field theory dynamics is captured by 
classical gravity as the dual string coupling is small and one moreover has a macroscopically 
large AdS$ spacetime. 

In the planar large N limit, the field theory dynamics truncates to the dynamics of 
single trace operators O = Tr (\&) where \& is a collection of the basic fields appearing in 
the Lagrangian. For M = 4 SYM one has the bosonic R-charge scalars transforming in the 
vector representation of the 577(4) R-symmetry, the gauge fields, and fermions; schematically 
$ = {X 1 , Dp, ip a }. The dynamics in the single trace sector in the planar limit is classical due 
to large N factorization, but the general structure is hard to derive even for the particularly 
symmetric case of M = 4 SYM. From a field theory viewpoint issues such a locality of 
this effective classical description are rather mysterious. However, once we pass to the 
holographic dual we encounter a manifestly local description in terms of two derivative 
gravity coupled to a bunch of matter fields in an asymptotically AdSs spacetime. In making 

19 The harmonic number function may be in fact be re-expressed in terms of the digamma function, or 
more simply as 

T '( x ) 

Harmomc(x) = 7 e + — — , 
T(x) 

where 7 e is Euler's constant; see [97] for a discussion of its properties. 

20 For four dimensional boundary field theories we can replace the S 5 by a Sasaki-Einstein manifold X5 
such as T ' or L pqr . In these cases one recovers jV = 1 quiver supcrconformal gauge theories on the 
boundary. 
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this statement we are assuming that we have performed a Kaluza-Klein reduction of the 
Type IIB supergravity fields over the compact S 5 leading to an infinite tower of massive 
fields coupled to the gravitational degrees of freedom. 

The general structure of this effective five- dimensional lagrangian is not only complicated 
but it also depends on the details of the internal space. Were one to replace the S 5 by a 
Sasaki-Einstein five manifold X 5 one would end up with a different effective description 
corresponding to a different field theory fixed point in four dimensions. However, there is a 
universal sub-sector of Type IIB supergravity which we can focus on - this is just the sector 
of solely gravitational dynamics in AdSs i.e., we set all the Kaluza-Klein harmonics of the 
graviton modes on S 5 and other matter degrees of freedom consistently to zero. We will 
restrict attention to this sub-sector which corresponds in the dual field theory to focussing 
on just the dynamics of the energy-momentum tensor. 

4.1 The universal sector: gravity in AdS^+i 

As discussed above we will concentrate on pure gravitational dynamics in asymptotically 
AdS spacetimes. This in particular allows us to work without loss of generality in arbitrary 
dimensions as the form the gravitational action is independent of the number of spacetime 
dimensions. Let us therefore consider starting with a string or M-theory background of 
the form AdS^+i xX where X is some compact internal manifold ensuring that one has a 
consistent string/M-theory vacuum. 21 The universal sector of this theory which we focus on 
is the dynamics of Einstein gravity with a negative cosmological constant, i.e., 



<Sbulk — 

16 



7TLr N J 



22 



With a particular choice of units (Ra<is = 1) Einstein's equations are given by 

Emn = Rmn — - GmnR — - — - Gmn = {A 2) 

Rmn + dG M N = 0, R=-d(d + l). 

Of course the equations (4.2) admit AdS^+i solutions, which correspond to the vacuum 
state of the dual field theory. Recall that global AdS^+i has as its boundary the Einstein 
static universe, Rx S d_1 . We are free to consider other choices of boundary manifolds Bd', for 
instance to discuss field theory on Minkowski space R d_1,1 we would focus on the Poincare 
patch of AdSd+i- Given a metric g on the boundary Bd we have the bulk geometry to zeroth 



21 We will be interested in d > 2. As discussed in [66, 69] there is no interesting hydro-dynamic limit for a 
1 + 1 dimensional CFT. A conserved traceless stress tensor is simply made up of left and right moving waves 
which propagate with no dissipation. 

22 We use upper case Latin indices {M, N, ■ ■ ■ } to denote bulk directions, while lower case Greek indices 
{/i, v, ■ ■ ■ } refer to field theory or boundary directions. Finally, we use lower case Latin indices {i,j, • • • } to 
denote the spatial directions in the boundary. 
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order in the Fefferman-Graham expansion given by 23 

ds 2 = — (dz 2 + dx» dx v ) . (4.3) 

We will refer to these spacetimes as asymptotically AdS^+i although this terminology is 
strictly speaking incorrect in a strict general relativistic sense. We will return to the issue 
of how to find solutions to (4.2) with this prescribed boundary metric, after discussing some 
issues related to fluid dynamics. 

We are interested in the description of the field theory in the canonical ensemble and 
in particular on situations where the field theory only attains local thermal equilibrium. 
There is usually an interesting phase structure for the field theories on non-trivial boundary 
manifolds Bd- Usually this arises from the fact that one can have dimensionless ratios of 
length scales arising from the non-trivial geometry of the background. The classic example 
is Bd = R x S d_1 where the low temperature phase is described to be the confined phase 
with O (1) free energy while the high temperature phase is the deconfined phase with O {N 2 ) 
free energy [98]. The former phase is dual to a thermal gas in AdS^+i while latter phase has 
a geometric description in terms of a Schwarzschild black hole in AdSd+i. 

To discuss hydrodynamics we need to be in the long- wavelength regime which is achieved 
only in the deconfined phase at high temperatures [65]. A simple way to see this is to 
note that in conformal field theories the phase structure is determined by the dimensionless 
ratio of length scales - if Bd has curvature scale R c and we are interested in the canonical 
ensemble at temperature T, then the phase structure depends on R C T. However, since the 
mean free path of the system £ m f p ~ 1/T from conformal invariance, we see that in order 
for the gradient expansion to be valid we require that R C T ^> 1. This can equivalently be 
interpreted as the statement that variations in the curvature of the background are small in 
units of the local temperature which means that we can approximate the boundary metric 
to be locally flat. This discussion is in perfect accord with the fact that for the CFT on 
Minkowski space the absence of any length scales in the background means that one has a 
trivial phase structure - the field theory is always deconfined on R :1 . 

Therefore to construct the dual of hydrodynamics on a curved boundary manifold Bd 
we can use as our starting point a flat boundary and systematically account for curvature 
terms as we proceed in the gradient expansion. In fact, the leading order viscous hydrody- 
namics is insensitive to the boundary curvature terms which only show up at second order 
in derivatives. 

23 In fact, if <?^„ is Ricci fiat then (4.3) satisfies the Einstein's equations (4.2). Similar constructions 
can be made for Einstein metrics on the boundary, albeit with a different warping function i.e., ds 2 = 
-7 [dz 2 + F(z) g^v dx 11 dx v ) for an appropriate choice of F(z). For instance F(z) = (1 + ^-) 2 for an Einstein 
manifold with negative curvature. 
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4.2 Preliminaries: Schwarzschild black holes in AdS^+i 

Let us consider the geometry dual to thermal field theory on Minkowski space, which is given 
by the planar Schwarzschild-AdS^+i black hole which in Schwarzschild type coordinates is 
given by 

dv 2 

ds 2 = -r 2 f(br) dt 2 + - — . + r 2 8^ dx l dx J , 

r 2 f{br) ^ 

While this is a one-parameter family of solutions labeled by the horizon size r + which sets 
the temperature of the black hole 

r = ^- < 4 - 5 > 

it is easy to generate a d parameter family of solutions by boosting the solution along the 
translationally invariant spatial directions x\ leading to a solution: 24 

dv 2 

ds 2 = ——— + r 2 {-f{br)u^u u + P^ u ) dx» dx u (4.6) 
r z j{br) 

with 



u 



u' 



(4.7) 



where the temperature T and velocities /3j are all constants with (3 2 = (3j ft , and P^ v = 
u^u u + r]^ is the projector onto spatial directions. These solutions are generated by a 
simple coordinate transformation from (4.4) and can be understood physically as follows. 
The isometry group of AdS^+i is SO(d,2). The Poincare algebra plus dilatations form a 
distinguished subalgebra of this group. The rotations SO(d) and translations R 1,d_1 that 
belong to this subalgebra annihilate the static black brane solution (4.4) in AdS^+i. This is 
manifest from the symmetries of the background (4.4). However, the remaining symmetry 
generators, which are the dilatations and boosts, act nontrivially on this brane, generating a 
d parameter set of black hole solutions. The parameters which characterize the bulk solution 
are precisely the basic hydrodynamical degrees of freedom, viz., temperature and the velocity 
of the black hole. 

The boosted black hole (4.6) is an asymptotically AdS^+i solution which has a holo- 
graphic stress tensor on the boundary. It is in fact not hard to see that the stress tensor 
for this solution is the ideal conformal fluid stress tensor (3.7) with a particular value of the 
normalization constant a = — ^ d+1) ■ This is not surprising as the solution is stationary 



The indices in the boundary are raised and lowered with the Minkowski metric i.e., = rj^ 
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and therefore corresponds to the global thermal equilibrium. To describe hydrodynamics we 
should perturb the system away from global equilibrium. Based on our discussion in §2 it is 
natural to consider a situation where the thermodynamic variables vary along the boundary 
directions. This can be simply achieved by promoting the parameters b, fa to functions 
of the boundary coordinates and furthermore letting the boundary metric vary to account 
for curvature couplings. As long as the variations are of large wavelength we can work in 
the effective field theory framework and construct a solution order by order in boundary 
derivative expansion. 

Roughly speaking, our construction may be regarded as the 'Chiral Lagrangian' for brane 
horizons. We have hitherto discussed how the dilatation and boost generators of the con- 
formal algebra act on the space of black hole solutions which is described by a point in d 
dimensional parameter space. Our construction effectively promotes these parameters to 
'Goldstone fields' (or perhaps more accurately collective coordinate fields) and determines 
the effective dynamics of these collective coordinate fields, order by order in the derivative 
expansion, but making no assumption about amplitudes. 25 

4.3 Regularity and choice of coordinate chart 

Whilst it is straightforward to promote the parameters b and fa appearing in (4.6) to func- 
tions depending on (t,x l ), there is a subtlety we need to consider. The issue is that the 
Schwarzschild coordinates used to write the metric in (4.6) are not regular on the future 
horizon. We would like to work with coordinates which are manifestly regular everywhere 
except for the singularity at the origin r = in (4.6), for we expect that the class of stress 
tensors that are fluid dynamical in nature should be special from a general relativistic view- 
point (we will shortly review why). In fact, we will show that the fluid dynamical stress 
tensors give rise to regular black hole spacetimes as their holographic duals. 

One can motivate the special nature of fluid dynamical form for the stress tensor by the 
following observation: consider the Fefferman-Graham form of the AdS^+i metric (4.3). It 
is well known that in order to find an asympotically AdS^+i spacetime with a prescribed 
boundary 3d, one has to give in addition to the metric g^ u on Ed another piece of data which 
corresponds to the boundary stress tensor. Armed with g^ v and T^ v one can construct a bulk 
solution order by order in a perturbative expansion in the Fefferman-Graham radial variable 
z. To leading order the solution is simply 

ds 2 = — (dz 2 + (g^ + az d T^) dx 1 * dx u ) . (4.8) 

z 

This scheme for constructing a bulk spacetime with prescribed boundary data is well devel- 
oped in the formalism of holographic renormalization [99] . However, this scheme is not likely 
to generically reproduce regular bulk spacetimes. To see why one just needs to a do a simple 

25 Another useful way to think about the hydrodynamic description in gravity language is to view it as 
the collective field theory of the lowest quasi-normal mode, the mode with vanishing frequency at zero 
momentum, of the black hole geometry. 
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degrees of freedom counting. A traceless, symmetric stress tensor on Bd has d ^ ^ — 1 degrees 
of freedom. But the dynamical equations of motion are simply the conservation equations 
(2.1) which are just d equations leading to a vastly underdetermined system when d > 2. 
However, a fluid dynamical stress tensor is a special class of conserved stress tensors for it 
is described by precisely d degrees of freedom, the temperature and velocity. 26 

In order to make regularity manifest, we will describe how to construct gravitational 
black hole solutions dual to arbitrary fluid flows using a coordinate chart that is regular on 
the future horizon. 27 We work with a set of generalized Gaussian null coordinates which 
are constructed with the aim of having the putative horizon located at some hypersurface 
r(x M ) = r#(x M ). 28 So as the starting point for our analysis we consider the boosted planar 
Schwarzschild-AdSrf + i black hole solution: 

ds 2 = -2 Ufl dx^dr - r 2 f(br) Ufl u v dx lL dx v + r 2 dx lL dx v , (4.9) 

where we have written the metric in ingoing Eddington-Finkelstein coordinates. We should 
note that it is possible to recast (4.9) in a Weyl covariant form when the boundary metric 
on Bd is curved - we have [69]: 

ds 2 = -2u^dx^ {dr + rA fl dx^) + r 2 (1 - f(br)) u^u y dx^dx v + r 2 g^dafdx" . (4.10) 

The main rationale behind switching to these Eddington-Finkelstein coordinates apart 
from making issues of regularity more transparent, is that they provide a clear physical pic- 
ture of the locally equilibrated fluid dynamical domains in the bulk geometry. The boundary 
domains where local thermal equilibrium is attained in fact extend along ingoing radial null 
geodesies into the bulk. So a given boundary domain corresponds to an entire tube of width 
set by the scale of variation in the boundary, see Fig. 1 for an illustration. In the Eddington- 
Finkelstein coordinates one just has to patch together these tubes to obtain a solution to 
Einstein's equations and moreover this patching can be done order by order in boundary 
derivatives, just as in fluid dynamics. We now proceed to outline a perturbation scheme 
which allows us to construct the desired gravity solution dual to fluid dynamics. 

As the starting point consider the metric (4.9) with the constant parameter b and the 
velocities replaced by slowly varying functions b(x IJ '), f3i(x fl ) of the boundary coordinates. 

ds 2 = -2 Ufl {x a )dx^dr-r 2 f{b{x a )r) Ufl {x a ) u v {x a ) dx» dx u + r 2 P^ v (x a ) dx» dx u . (4.11) 

Generically, such a metric (we will denote it by (b(x^) , Pi(x fJj )) is not a solution to Ein- 
stein's equations. Nevertheless it has two attractive features. Firstly, away from r = 0, 
this deformed metric is everywhere non-singular. This pleasant feature is tied to our use of 
Eddington-Finkelstein coordinates. Secondly, if all derivatives of the parameters b{x^) and 

26 Similar arguments can be given if we want to consider fluids which carry conserved global charges. 
27 For a discussion on the Fefferman-Graham coordinates and regularity see [100]. 

28 We will determine explicitly the location of the horizon after sketching the construction of the solution. 
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Fig. 1: Penrose diagram of the uniform black brane and the causal structure of the spacetimes dual to 
fluid mechanics illustrating the tube structure. The dashed line in the second figure denotes 
the future event horizon, while the shaded tube indicates the region of spacetime over which 
the solution is well approximated by a tube of the uniform black brane. 

Pi(x^) are small, is "tubewise" well approximated by a boosted black brane. Conse- 
quently, for slowly varying functions it might seem intuitively plausible that 
(4.11) is a good approximation to a true solution of Einstein's equations with a regular event 
horizon. In [4] this intuition is shown to be correct, provided the functions 6(0^) and j3i{x^) 
obey a set of equations of motion, which turn out simply to be the equations of boundary 
fluid dynamics. 

Einstein's equations, when evaluated on the metric G^°\ yield terms which involve deriva- 
tives of the temperature and velocity fields in the boundary directions (i.e., (xi,v) = x^) 
which we can organise order by order in a gradient expansion. Note that since is an 
exact solution to Einstein's equations when these fields are constants, terms with no deriva- 
tives are absent from this expansion. It is then possible to show that field theory derivatives 
of either \nb(x fl ) or f3i(x^) always appear together with a factor of b. As a result, the con- 
tribution of n derivative terms to the Einstein's equations is suppressed (relative to terms 
with no derivatives) by a factor of (b/L) n ~ 1/(T L) n . Here L is the length scale of varia- 
tions of the temperature and velocity fields in the neighbourhood of a particular point, and 
T is the temperature at that point. Therefore, provided LT ^> 1, it is sensible to solve 
Einstein's equations perturbatively in the number of field theory derivatives. 29 Essentially 
we are requiring that 



29 Note that the variation in the radial direction, r, is never slow. Although we work order by order in 
the field theory derivatives, we will always solve all differential equations in the r direction exactly. This 
should be contrasted with the holographic renormalization group which is a perturbativc expansion in the 
FcfFcrman- Graham radial coordinate [99]. 




(4.12) 
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where we introduce a book-keeping parameter e which keeps track of derivatives in the 
boundary directions. It is useful to regard b and j3i as functions of the rescaled field theory 
coordinates ex^ where e is a formal parameter which we eventually set to unity. This way 
every derivative of $ or b produces a power of e, consequently powers of e count the number 
of derivatives. 



4.4 The perturbative expansion in gravity 

We now describe a procedure to solve Einstein's equations in a power series in e. Consider 
the metric 30 

G = b) +eGW(/3 i , b) + e 2 G^(A, b) + O {e 3 ) , (4.13) 

where G^ is the metric (4.11) and G^ etc., are correction metrics that are yet to be 
determined. We will see that perturbative solutions to the gravitational equations exist only 
when the velocity and temperature fields obey certain equations of motion. These equations 
are corrected order by order in the e expansion; this forces us to correct the velocity and 
temperature fields themselves, order by order in this expansion. Consequently we set 

A = A (0) + e 0® + O {e 2 ) , b = &<°> + e&W + O {e 2 ) , (4.14) 

where and b^ are all functions of e x M . 

In order to proceed with the calculation, it will be useful to fix a gauge. In [4] it was 
convenient to work with the 'background field' gauge 

G rr = , G rfl oc , Tr ((G (0) ) _1 G (n) ) =0 Vn>0. (4.15) 

One could equivalently work with a slightly different gauge choice [69]: 

G rr = , G rii = n M . (4.16) 

With this gauge choice it transpires that curves of x M = constant are in fact affinely parame- 
terized null geodesies in the resulting spacetime, with the radial coordinate r being the affine 
parameter. With the former gauge choice (4.15) of course = constant are null geodesies; 
however, r is not an affine parameter for this geodesic congruence. 

With this picture in place one can plug in the ansatz (4.13) and (4.14) into Einstein's 
equations (4.2) and expand them order by order in e. Let us imagine that we have solved 
the perturbation theory to the (n — l) th order, i.e., we have determined G^ for m < n — 1, 
and have determined the functions /3 t - and b^ for m < n — 2. Plugging the expansion 
(4.13) into Einstein's equations, and extracting the coefficient of e n , we obtain an equation 
of the form 

G (0) ($° } , & (0) )1 G (n) (^) = Sn . (4.17) 



30 For convenience of notation we are dropping the spacetime indices in G^ n > . We also suppress the 
dependence of b and on x M . 
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Here HI is a linear differential operator of second order in the variable r alone. As is 
already of order e n , and since every boundary derivative appears with an additional power 
of e, H is an ultralocal operator in the field theory directions. In fact not only is EI is a 
differential operator only in the variable r independent of but also its precise form at the 
point x M depends only on the values of fif^ and at x M , and not on the derivatives of these 
functions at that point. Furthermore, the operator H is independent of n; we have the same 
homogeneous operator at every order in perturbation theory. This makes the perturbation 
expansion in e ultra-local in the boundary directions; we can solve the equations point by 
point on the boundary! 

The source term s n however is different at different orders in perturbation theory. It 
is a local expression of n th order in boundary derivatives of (3^ and b^°\ as well as of 
(n — k) th order in for all k < n — 1. Note that (3^ and do not enter the n th 

order equations as constant (derivative free) shifts of velocities and temperatures solve the 
Einstein's equations. 

The gravitational equation (4.17) form a set of ( d+1 *)( d+2 *) equations. It is useful to split 
these into two classes of equations: (i) a class that determines the metric data we need, 
comprising of f^±Ii equations which we view as dynamical equations and (ii) a second set of 
d equations which are essentially constraint equations. 



Constraint equations: We will refer to those of the Einstein's equations that are of 
first order in r derivatives as constraint equations. These are obtained by contracting the 
equations with the one-form normal to the boundary 

E$ = E MN £ N (4.18) 

where for our considerations £/v = dr. Of these equations, those with legs along the boundary 
direction are simply the equations of boundary energy momentum conservation: 

V^ ( £i) = • (4.19) 

Here T^_ x s is the boundary stress tensor dual the solution expanded up to O (e n_1 ) and is a 
local function of the temperature and velocity fields involving no more than n — 1 derivatives. 
Furthermore, it is conformally covariant and consequently it is a fluid dynamical stress tensor 
with n — 1 derivatives. 

These constraint equations can be used to determine M n_1 ) and /?f n ^; this is essentially 
solving the fluid dynamics equations at order O (e n ) in the gradient expansion assuming that 
the solutions at preceding orders are known. There is a non-uniqueness in these solutions 
given by the zero modes obtained by linearizing the equations of stress energy conservation 
at zeroth order. These can be absorbed into a redefinition of the f3^ , b^ , and do not 
correspond to a physical non-uniqueness. 
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Dynamical equations: The remaining constraint E rr and the dynamical Einstein's equa- 
tions E^ v can then be used to solve for the unknown function g^ n \ By exploiting the un- 
derlying symmetries of the zeroth order solution, specifically the rotational symmetry in the 
spatial sections on the boundary, SO(d— 1), it is possible to decouple the system of equations 
into a set of first order differential operators. Having performed this diagonalization of the 
system of equations one has a formal solution of the form: 

G^ n ' = particular (s n ) + homogeneous (H) (4.20) 

To determine the solution uniquely we need to prescribe boundary conditions: we impose that 
our solution is normalizable so that the spacetime is asymptotically AdS^+i and also demand 
regularity at all r / 0. In particular, the solution should be regular at the hypersurface 
br = 1. It has been shown in [4] that for an arbitrary non-singular and appropriately 
normalizable source s n encountered in perturbation theory, it is always possible to choose 
these constants to ensure that is appropriately normalizable at r = oo and non-singular 
at all nonzero r. Furthermore, if the solution at order n — 1 is non-singular at all nonzero 
r, it is guaranteed to produce a non-singular source at all nonzero r. Consequently, the 
non-singularlity of s n follows inductively. 31 



4.5 Details of the long-wavelength perturbation expansion 

We have described how to perturbatively construct solutions to (4.2) by working order by 
order in boundary derivatives. We now briefly illustrate how to carry this out in practice to 
first order in derivatives and refer the reader to the original references [4, 69] for the detailed 
derivation of the results given herein. 

Consider the zeroth order metric G^ given in (4.11). If we want to work to first order 
in boundary derivatives, we can pick a point on the boundary x M = Xq, which by exploiting 
the Killing symmetries of the background can be chosen to be the origin. At Xq we can use 
the local scaling symmetry to set = 1 and pass to a local inertial frame so that = 0. 
Expanding (4.11) to O (e) at this point we find 

dsj j = 2 dv dr — r 2 f(r) dv 2 + r 2 dx^ dx % 

- 2 5pf ] dx 1 dr-2 5pf ] r 2 (1 - f(r)) dx 1 dv - d ^ dv 2 , 

where we have introduced a short-hand 5(3^ = x M d^fif^ and 5b^ = x M d^b^ which are the 
leading terms in the Taylor expansion of the velocity and temperature fields at Xq (taken to 
be the origin). 

We now need to pick an ansatz for the metric correction at O (e), G^\ which we wish 
to determine. As mentioned earlier it is useful to exploit the SO(d — 1) spatial rotation 

31 There is a slight subtlety which needs to be borne in mind here: the requirements above do not completely 
fix G(") since the differential operator H has some zero modes. These can however be fixed by appropriately 
absorbing the zero modes into redefinition of the local temperature and velocity fields. 
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symmetry at Xq to decompose modes into various representations of this symmetry. Modes 
of G^ 1 ' transforming under different representations decouple from each other by symmetry. 
We have the following decomposition into SO(d — 1) irreps: 



scalars: ? ? ^ 



(i) 
a 5 



vectors: G vi 
tensors: G$ 



(i) (4.22) 



We work sector by sector and solve the constraint and the dynamical Einstein's equations. 
In the scalar sector, we find that the constraint equations imply that 

d4 0) = d v b® , (4.23) 



d- 

while in the vector sector we have 

= d v pl 0) . (4.24) 

These two equations are simply the equations of energy momentum conservation (2.1) at the 
point Xq. The remaining dynamical equations are to be solved for the functions appearing 
in G^ - we refer the reader to [4, 69] for the specifics and just record here the form of the 
differential operator we obtain in various sectors: 32 

d ( 1 d ,„\ 
vector : M^O = —(—- — O ) 

dr \r a 1 dr J 

tensor : I d(d+1) = 4" ( r<m f( r ) 4 O) (4.25) 
2 dr \ dr J 

which as advertised earlier are simple differential operators in the radial variable alone and 
clearly can be inverted to find the function O once the source s n is specified. 

The calculation can in principle be carried out to any desired order in the e expansion. 
As we have discussed earlier the form of the differential operator (4.25) remains invariant in 
the course of the perturbation expansion. What one needs to compute at any given order 
are the source terms s n . In addition one has to always ensure that the lower order stress 
tensor conservation equations are satisfied. For instance in order for the source terms which 
appear in the determination of G^ at second order to be ultra-local at our chosen boundary 
point Xq we have to ensure that the first order fluid equations of motion are satisfied. This 
is encapsulated in our discussion of the constraint equations in the bulk (4.19). 



32 We have indexed the operator H by the representation label of the SO(d — 1) rotational symmetry. The 
scalar sector involves some mixing between different fields and is slightly more involved, see [4, 69] for details. 
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5 Gravitational analysis: Metrics dual to fluids 



We have thus far discussed how to solve to Einstein's equations order by order in boundary 
derivatives. We now present the result for the general fluid dynamical metric up to second 
order in boundary derivatives and then describe how one extracts the stress tensor quoted in 
§3.4. Subsequently, we will analyze the physical properties of these solutions and argue that 
they are regular black hole solutions with an event horizon. Following that we will discuss 
how one can use the black hole nature of the solution to understand aspects of the fluid 
dynamics such as the entropy current. 

5.1 The gravitational dual to non-linear viscous fluid 

We have described how one can consistently solve for the bulk metric in §4 - in particular, 
in §4.5 we have given a sketch of how the perturbation scheme works. It turns out that the 
bulk metric resulting from the explicit computation can schematically be cast into the form: 

ds 2 = G MN dX M dX N = -2S(r,x)u ll (x)dx^dr + x^(r,x)dx'"dx u . (5.1) 

This was the form of the metric originally derived in [4], where the functions S(r,x) and 
x(f, x) are explicit functions of the radial coordinate r, whilst being given in terms of a 
gradient expansion in the boundary directions x^. Note that the bulk metric is actually 
given in the gauge (4.15). The expressions for the functions S(r,x) and x( r , x ) are rather 
cumbersome and we will not record them here, but rather refer the interested reader to the 
original source [4, 66]. Instead we will record the explicit metric below making manifest 
the Weyl covariant structures, following [69] which makes for a more compact presentation. 
Before we do that however, it is useful to understand the geometry in the somewhat simpler 
form of the metric given above, since it captures much of the essential physical features. 

First of all, it is useful to realize that in the metric (5.1) lines of constant x M are radially 
ingoing null geodesies. If the function S(r, x M ) is gauge fixed to unity then in fact the radial 
coordinate r is actually an afflne parameter along this null congruence. One can visualize this 
as follows: consider a null geodesic congruence i.e., a family of null geodesies with exactly 
one geodesic passing through each point, in some region of an arbitrary spacetime. Let E be 
a hypersurface that intersects each geodesic once. Let x M be coordinates on E. Now ascribe 
coordinates (p, x M ) to the point at an affine parameter distance p from E, along the geodesic 
through the point on E with coordinates x M . Hence the geodesies in the congruence are lines 
of constant x^. In this chart, this metric takes the form 

ds 2 = -2 u^x) dp dx^ 1 + x^ip, x) dx^ dx v ', (5.2) 

where the geodesic equation implies that is independent of p. It is convenient to generalize 
slightly to allow for non-affine parametrization of the geodesies: let r be a parameter related 
to p by dp/dr = S(r, x). Then, in coordinates (r, x), the metric takes the form given in (5.1). 
Note that E could be spacelike, timelike, or null. We shall take E to be timelike. 
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The metric (5.1) has determinant —S 2 x' Iu u fl Updetx, where x^ u is the inverse of Xiwi 
where the indices are raised with the boundary metric; for (5.1) the induced metric on the 
boundary Bd is simply 

= lim -^x^uir.x) . 

r^oo r 

Hence the metric and its inverse will be smooth if S, u^ and Xtiv are smooth, with 5^0, x^v 
invertible, and x^ v u \i timelike. These conditions are satisfied on, and outside, the horizons 
of the solutions that we shall discuss below. Finally, note that the inverse metric, 33 G MN , 
can be determined easily to be 

S X aP up QQ& _ S 2 u, u 5 ( X aP X ,S ~ X ai X PS ) 



G 



S 2 Un Up X 



[IV 



QTC 



S 2 Un U v X 



[IV 



-S 2 Ujj_ u v x 



(5.3) 



using G M k G 



KN 



ft N 
°M ■ 



Weyl covariant form of fluid metric: As remarked above the coordinates in which we 
present the fluid metric (5.1) do not make explicit the Weyl covariant structures. An alternate 
form of the metric was written down in general for asymptotically AdS^+i spacetimes in [69] 
which takes the form: 



ds 2 = G M n dX M dX N = -2 Ufl (x) dx" (dr + Q3 v (r, x) dx u ) 



(5 tiV (r,x)dx"dx u , (5.4) 



where the fields 5J M and (5^ are functions of r and x^ which admit an expansion in the 
boundary derivatives. In the parameterization used in [69] one finds the metric functions are 
given up to second order in derivatives as: 
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where we are using the tensors defined earlier in (3.23) and also introduce the Schouten 
tensor (a particular combination of curvature tensors) which is defined in (A. 6). The 
tensor (5 M1/ is clearly transverse, since it is built out of operators that are orthogonal to the 
velocity, and it can be inverted via the relation ((S" 1 )' 1 " & au = P%. We also note that the 
coordinate chart used to write (5.1) is consistent with the bulk gauge choice (4.16). For 
future reference we record the induced metric on the boundary in these coordinates 
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Note that the 'inverse d-metric' x 111 ' is defined via Xfiv X vp = S.f. 
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which is crucially used to raise and lower the boundary indices (lowercase Greek indices). 
The various functions appearing in the metric are given in terms of definite integrals 
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The asymptotic behaviour of these functions 0i(y) and K>i(y) which are relevant for the stress 
tensor computation can be found in [69]. 

5.2 The boundary stress tensor 

Once the bulk black hole solution is determined it is straightforward to use the holographic 
prescription of [95, 96] to compute the boundary stress tensor. To perform the computation 
we regulate the asymptotically AdS^+i spacetime at some cut-off hypersurface r = A c and 
consider the induced metric on this surface, which up to a scale factor involving A c is our 
boundary metric g^ v . The holographic stress tensor is given in terms of the extrinsic curva- 
ture and metric data of this cut-off hypersurface. Denoting the unit outward normal to 
the surface by n M we have 

^ = S W V\ (5.8) 
For example for asymptotically AdSs spacetimes the prescription of [96] gives 
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-Kg^-id- 1) g*™ - -^—^ - -Rg^ 



(5.9) 



where K^ u is the extrinsic curvature of the boundary. Implementing this procedure for the 
metric (5.4) and utilizing the asymptotic form of the functions given in (5.7) we recover the 
stress tensor quoted in (3.23) with the precise transport coefficients (3.25). 

5.3 Event horizons 

Having understood the geometric aspects of the coordinate chart employed for metrics dual 
to arbitrary fluid flows in the boundary, we are now in a position to address our assertion 



36 



that these metrics are in fact black hole spacetimes with regular event horizon. In general 
determining the event horizon of a black hole spacetime requires that we know the entire 
future evolution of the spacetime geometry. This is basically due to the definition of an event 
horizon, which is teleological in nature. 

Formally one defines the event horizon of a given spacetime as follows: the future event 
horizon 7i + is the boundary of past set of future null infinity. This is the mathematical 
way to capture the physical statement that the spacetime events inside the event horizon of 
the black hole cannot communicate to the asymptotic region. It is important to note that 
the future null infinity J" + corresponds to where the future directed null geodesies in the 
spacetime end and is in fact timelike for asymptotically AdS spacetimes. Since 7i + is the 
boundary of a causal set, it is a null surface which is in particular generated by null geodesies 
in the spacetime. One can determine the event horizon by shooting null geodesies from the 
interior of the spacetime and checking whether they make it out to J^ + . By fine tuning the 
initial conditions of the geodesic one can zero onto 7i + . Alternately, if one knows the late 
time generators of the event horizon, then one can evolve the geodesic equation backwards 
with these generators as the future boundary condition. 

For the metrics dual to boundary fluid dynamics, we will make an important assumption: 
since fluid dynamical evolution tends generically to smooth out inhomogeneities, we will 
assume that the late time solution is one corresponding to global equilibrium wherein the 
fluid settles down. Then at late times we have a clear idea of where the event horizon is 
located - its location in the radial direction in the asymptotically AdS^+i bulk spacetime is 
specified by the value of the local temperature. Armed with this data one can in principle 
evolve the null generators of the horizon backwards to construct the entire surface 7i + . 34 
However, it will turn out that in the long- wavelength approximation one can as usual work 
order by order in the boundary derivative expansion and in fact determine the location of 
the horizon. 

For the spacetime (5.1), let us suppose that the null hypersurface corresponding to 7i + 
that we are after is given by the equation 

S<n(r,x) = , with Sn(r,x) = r — rn(x) . (5.10) 

As we are working in a derivative expansion and assuming that the dissipative fluid dynamics 
drives us towards global equilibrium with a local temperature T(x) = ^rf^y we take 

^(x) = ^y + ^£ fc r (fc )(x). (5.11) 

where r^(x) denote the corrections away from the late time hypersurface in the spacetime. 

34 As has been recently discussed in the context of dynamical fluid flows [63] there are examples such as the 
the Bjorken flow and the conformal soliton flow where despite the long wavelength approximation required 
for fluid dynamics to be valid being upheld, one nevertheless has non-trivial late time boundary conditions 
which change the nature of the event horizon. 
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By demanding that the hypersurface defined by (5.10) be null i.e., 

G AB d A S n d B S n = (5.12) 

we can solve for the functions r^(x) order by order in the e expansion. In addition to the 
location of the horizon we will also be interested in the normal vector to the event horizon 
by £ A : by definition, 

i A = G AB d B S H (r,x) (5.13) 

which also has an e expansion. Using the explicit form of the metric (5.1) or (5.4) one 
can easily write down the equation determining the location of the event horizon (5.12) to 
arbitrary order in e explicitly. For instance in the simpler situation of the non-Weyl covariant 
metric (5.1) we have the equation 

= ^ {l-2eS X aP up d a r H - e 2 S 2 ( X a ? X 1 ' ~ X° 7 X P8 ) *h u s d a r H d P r H ) , 

(5.14) 

which is actually an algebraic equation for the functions rn^ at any given order in e and 
hence easy to solve. For the Weyl covariant metric (5.4), the horizon location is determined 
by the equation [69]: 

(0- 1 ) f k (J k„ = 2u' , ^ 1 (5.15) 

with 

= d^n + ^3^(r H (x),x) , (5.16) 

where we have left the powers of e implicit. In terms of k m one can determine the components 
of the normal vector to the horizon hypersurface as 

£ A d A = n" (d„ + dpr H d r ) , with = - {p-^ k v . (5.17) 

Rather than illustrate the computation of the event horizon for the metrics (5.4) or (5.4) 
we will just give the resulting answer of the computation, referring the reader to the original 
derivation presented in [66, 69]. However, before doing so we would like to present a toy 
model computation which illustrates the key features of the computation. 

Toy model for event horizon detection: The general situation we are interested in has 
dependence on all the boundary directions x^ - a simpler situation to consider would be 
to look at a case where we have dependence only on one variable, say time. Furthermore, 
there is nothing really special about asymptotically AdS spacetimes for the detection of the 
event horizon. One could just as readily have made similar arguments for asymptotically 
flat spacetimes. As a result we will focus on a simple example of the Vaidya spacetime which 
describes a spherically symmetric black hole with infalling null matter. As we are interested 
in just the geometric properties of the spacetime, we will just focus on the metric and not 
worry about solving Einstein's equations. 
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Consider the Vaidya spacetime, whose metric is given as: 

/ 2m(v)\ 

ds 2 = - 1 — ) dv 2 + 2dvdr + r 2 dtt 2 . (5.11 



r 

We want to determine the location of the horizon, which by spherical symmetry has to lie 
on the locus r = r(v). The normal one-form to this surface is clearly just n = dr — rdv. 
Demanding that this be null gives us a differential equation for the null surface 

r(v) = 2 m(v ) + 2 r(v) r(v ) , (5.19) 

which is the analog of (5.14) for the Vaidya metric (5.18). The equation we have at hand 
is a first order ODE for the function r(v), solving which we would determine the location 
of the horizon non-locally in terms of m(v). This would of course be the case for a generic 
function m(v). 

However, to make contact with our earlier discussion, suppose we assume that the mass 
function m(v) is slowly varying and moreover that it asymptotes to a constant for late times, 
i.e., 

rh(v) = O (e) ,mm = 0(e 2 ) etc., and lim m(v ) = mo (5.20) 

v— too 

then we can solve (5.19) in a derivative expansion. Consider the ansatz, 

OO 

r(v) = 2m(v) + ^ e k r (k) (v) (5.21) 

k=l 

which leads to the solution 

?"(!)= 8mm, r ( 2 ) = 64m m 2 + 32 mm 2 , etc.. (5.22) 

Thus, by invoking the slow variation ansatz we obtain a local expression for the location of 
the horizon in a derivative expansion. 

Horizon location for fluid metrics: An analogous analysis can be carried out for the 
the equation (5.15) - we quote here just the final result (once again dropping e for brevity) 

r H (x) = 77^ + b(x) (Ni a aP a af3 + N 2 uj a/3 uo a(3 + K 3 K) (5.23) 

0{X) 

with 

2(d 2 + d-4) 2o 2 (l) 



d 2 (d-l)(d-2) d(d-l) 



d(d-l){d-2) 

TZ here is the Weyl covariant boundary Ricci scalar, which is defined in (A. 4). This establishes 
the claim made in the §1 that spacetimes dual to hydrodynamical evolution in the boundary 
field theory are regular black hole spacetimes. See Fig. 2 for an illustration of the event 
horizon in the spacetime. 
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Fig. 2: The event horizon r = r#(x M ) sketched as a function of the time t and one of the spatial 
coordinates x (the other two spatial coordinates are suppressed). 

5.4 Boundary entropy current 

As we discussed in §2.2 the flow of a dissipative fluid is characterized by entropy production. 
Given that we have constructed a gravitational dual which is a black hole it is natural to 
ask whether one can associate the entropy of the fluid with a geometric feature in the bulk 
spacetime. There is a natural object which one is tempted to use for this purpose, which is 
the area of the event horizon in the bulk geometry. 

In general when we consider deviations from equilibrium it is a-priori not clear that 
there should exist an unambiguous notion of entropy. In fact, since all we require of the 
entropy current is to satisfy the second law (2.20), we can use any local function having 
positive divergence to characterize the irreversibility of the fluid dynamical flow. The only 
constraint we would demand on this putative Boltzmann H-function is that it agrees with the 
thermodynamic notion of entropy in global equilibrium. For a stationary black hole, which 
corresponds to global thermal equilibrium, the area of the event horizon does indeed capture 
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the entropy of the dual field theory. It seems therefore natural to associate the entropy of 
the field theory with the area of the event horizon. This actually turns out to be a bit more 
subtle as discussed recently in [63]- we will return to this point after providing a sketch of 
the argument given in [66] for constructing an entropy current dual to fluid flows using event 
horizon data. 35 

Entropy current from geometry: We have established that the spacetimes dual to 
boundary fluid dynamics are in fact black hole geometries with a regular future event horizon 
7i + . Given a black hole spacetime one can associate an entropy with the area of the event 
horizon, which is the usual Bekenstein-Hawking entropy of the back hole. However, the field 
theory lives on the boundary of the spacetime and we would like to define an entropy current 
directly in the field theory. To this end, consider spatial sections of the event horizon, which 
are co-dimension two surfaces in the spacetime, which we label as 7i+. We are working 
in a coordinate chart where the coordinates a 1 for i = {1, • • • d — 1} define a chart on the 
spatial section and we use as affine parameter, the boundary coordinate v, to propagate these 
surfaces forward along the horizon generator £ A . On the surface 7i+ it is natural to define 
an area (d — l)-form, whose integral gives the area of the spatial section. From the area 
theorem for black holes, it follows that the area of the spatial sections will be non-decreasing 
with v . 36 

Once we have a geometric object such as the area-form on the horizon, we can pull it 
back to the boundary and have a candidate entropy function, for the pull-back too will 
have non-negative divergence. The only issue is the precise manner in which we implement 
the pull-back. It turns out to be natural to pull-back the area form on the horizon using 
radially ingoing null geodesies, which provide an isomorphism between the spatial sections 
on the boundary and the corresponding 7i+ on 7i + . This procedure was described in detail 
in [66] and was used there to construct a boundary entropy current explicitly. One can give 
a compact formula for the boundary entropy current based on this discussion as follows. Let 
us suppose that the induced metric on the spatial slices of the horizon TC^ is given by rjjj 
and as usual the boundary metric is g^ v . Under the split of the future horizons into spatial 
sections 7i+ we can also split the vector defined in (5.17) which essentially gives us the 
projection of the null generator of the event horizon into the boundary directions. Armed 
with this set of geometric data we define an entropy current 



which implements the pull-back procedure from the event horizon to the boundary. 

35 The rationale behind postponing the discussion of the subtleties is that the event horizon serves to 
illustrate the general idea of defining a boundary entropy current; it is trivial to generalize the construction 
to other quasi-local horizons. 

36 We assume that the null energy condition is satisfied. This is clearly true of the lagrangian (4.1), but 
we will demand the same to be true for the general discussion (6.1). 
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Boundary entropy current: The entropy current in d dimensions has to be a Weyl 
covariant vector of weight d. This follows from the fact that the entropy density scales like the 
inverse spatial volume, the total entropy being dimensionless, and the scaling of the velocity 
field given by(3.2). In [84] the general constraints on entropy currents were analyzed in the 
Weyl covariant formalism and a particular entropy current for four dimensional CFTs was 
constructed. This analysis was slightly generalized in [66] to argue for a five parameter family 
of fluid dynamical entropy currents in four dimensional CFTs. A similar analysis was carried 
out in arbitrary dimensions in [69] , who found a four parameter family of hydrodynamical 
entropy currents. The extra parameter which appears in d = 4 is related to the fact that one 
can have pseudo- vector contributions to the entropy current involving £ M defined in (2.16). 
In general d dimensions the entropy current takes the form: 37 

r s =su» + sb 2 u» (A 1 a aP a al3 + A 2 u aP ^ + A 3 Tl) 
+ s b 2 (B 1 V x a^ + B 2 2W A ) + • • • 

where s is the entropy density and ^1,2,3, -81,2 are arbitrary numerical coefficients. Requiring 
positivity of the divergence one finds a single linear relation between two of the coefficients: 

B 1 + 2A 3 = 0. (5.27) 

The gravitational analysis described above leads to a specific entropy current, with the 
parameters fixed by the geometric data encoded in the metric (5.4). The essential data is 
captured the entropy density, which in general is given by the Bekenstein-Hawking formula 
in terms of the area of the event horizon 38 

11 ' d (5.28) 
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The coefficients appearing in (5.26) are given in terms of the numerical (dimension) depen- 
dent constants satisfying (5.27) 

A 1 = ^(d + 2)- fI fl2 (l) + HtJ 2 (l) 
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Furthermore, it is possible to show that one can write the divergence of this gravitational 
entropy current as 



v p j s = ~y y 7 ^ + 2 
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[1 + A 1 d)-r w 



u a V a a^ +■•• , (5.30) 



37 The contribution from the pseudo- vector in four dimensions takes the form C\ sbi^ + C% s b 2 u x T> . 
The positivity condition on the divergence of the entropy current demands C\ + C2 = 0. 

38 As discussed above in the situations where the temporal variation is suitably slow the ateleological 
behaviour of the event horizon doesn't play an important role. Furthermore, the various quasi-local horizons 
are 'sufficiently near' the actual event horizon which makes it possible to use the area of the event horizon 
for our purposes. 
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which is accurate up to third order in the derivative expansion and clearly satisfies the 
requirement of non-negative divergence to that order. We have recast this expression purely 
in terms of fluid dynamical variables and it generalizes the result for viscous relativistic fluids 
(2.21). 

Event horizon vs. quasi-local horizons: We now turn briefly to the subtlety mentioned 
at the beginning of this sub-section with regard to using the area of the event horizon for 
capturing the entropy of the dual field theory in situations out of global equilibrium. A key 
feature of the event horizon is its teleological nature, i.e., the fact that one needs to know 
the entire future evolution of the spacetime in order to determine its location. This aspect of 
event horizons generically implies that one typically may encounter a horizon in the bulk even 
before we perturb the system - the event horizon grows in anticipation of the perturbation 
one is about to impart. Associating an entropy current with the event horizon as a result 
leads to a non-local and acausal definition of entropy as has been recently noted in [101] 
(cf., [63]). From an underlying statistical description, we would like however to define an 
entropy current a la Boltzmann, a local H-function which is defined purely based on the local 
dynamics of the fluid, with no recourse to late time boundary conditions. Moreover, there 
is a simple hydrodynamic flow on R" -1 ' 1 , the conformal soliton flow [53], where it has been 
shown that the event horizon area does not capture the entropy of the dual field theory [63]. 
Based on the latter analysis, it appears that in in certain dynamical situations one should 
use the area of apparent horizons (more precisely dynamical horizons) , to define the entropy 
current. 39 However, this statement glosses over certain subtleties involved in defining such 
quasi- local horizons 40 - for an account of the issues involved we refer the reader to [63]. 

6 Generalizations of the fluid-gravity correspondence 

In §4 and §5 we have discussed how one can construct gravitational solutions dual to fluid 
dynamics concentrating just on energy- momentum flow. One of the advantages of focussing 
on the stress tensor T^ u is that in the holographic context one deals purely with gravitational 
dynamics in AdS^+i. This allows for an universal description; the Lagrangian (4.1) is a 
consistent truncation of the bulk string theory in the AdS/CFT context (all matter field 
interact at best quadratically with the gravitational field). Nevertheless, one can incorporate 
other fields into the fluid-gravity correspondence and by now there is a large body of literature 
exploring these issues as discussed at the end of §1. We will describe the key features of 
these investigations using an abstract model and suggest some future generalizations. 

39 Another piece of evidence in favour of quasi-local horizons comes from the analysis of holographic 
entanglement entropy (for specified regions on the boundary) in time dependent states of the dual field 
theory, see the discussion in [102]. 

As explained in [63] we use the term apparent horizon to denote a co-dimension one surface in the 
spacetime, which in conventional general relativity literature would correspond to a marginally outer trapped 
tube. 
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6.1 Fluid-gravity and the inclusion of matter 

Consider a situation where we have a gravitational action which is described by the Einstein- 
Hilbert action with negative cosmological constant coupled to some matter fields which we 
collectively denote as To be precise, 

I d d+l x V^G (R-2A) + £ matter Gmn) (6.1) 
G N J 

where we have included metric contributions into the matter Lagrangian explicitly to allow 
for situations where we consider higher derivative theories of gravity (say for example the 
Gauss-Bonnet term in d > 4 as discussed in [79]). 

We will assume that the action (6.1) admits stationary asymptotically AdS^+i black 
hole solutions 41 which correspond to global thermal equilibrium, which by picking suitable 
coordinates we can write as 42 

ds 2 = -2 S(r, Q) Up dx 11 dr - r 2 V(r, Q) u u dx^ 1 dx v + r 2 cfe M dx v 

(6.2) 

where we have once again resorted to the ingoing Eddington-Finkelstein coordinates. The 
coordinate x^ as before span the boundary directions with \x = {0, • • • , d — 1}. Furthermore, 
we assume that we have translational invariance in the spatial directions which we have 
exploited to boost the black hole solution so as to allow for the velocity field. This of course 
requires that the matter field \& supporting the black hole to respect the symmetries of the 
metric. We have captured all the parameters describing the black hole into a single variable 
Q, which enters into the metric and the matter field and corresponds to all the physical 
charges we wish to ascribe to the geometry, such as temperature T, Maxwell charges qj 
(or equivalently chemical potentials fii), etc.. Depending on the boundary conditions on 
the matter field ^> we would find ourselves working in the canonical or the grand canonical 
ensemble. 43 It is possible to relax some of the assumptions described above at the expense 
of complicating the discussion. 

To obtain a geometry dual to fluid dynamics we want to promote the parameters of the 
background to fields as in conventional collective coordinate quantization. In order to do 
so we need to identify the full set of zero-modes. It is worthwhile to pause to take stock of 
some of the examples that have been discussed in the literature to orient ourselves: 



41 The general scheme we describe below is not restricted to asymptotically AdS spacetimes. The only 
reason to focus on these cases is the gauge-gravity duality which allows us to relate the bulk gravity dynamics 
to that of a boundary field theory and in particular allows us via the holographic renormalization scheme to 
extract a boundary stress tensor. One can construct inhomogeneous dynamical black holes using the slow 
variation ansatz for any desired asymptopia (including asymptotically flat spacetimes). We will discuss an 
example with different asymptopia in §6.2. 

42 In the following discussion we eschew the use of the Weyl covariant form since our disucssion would also 
be applicable to non-conformal fluids. 

43 Recall that for charged black holes which corresponds to the situation where '5 is a bulk Maxwell field, 
the grand canonical ensemble corresponds to non-zero value for the scalar potential Aq which is the boundary 
chemical potential fx. 
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1. In the forced fluid dynamics discussed in [70], ^> is just the dilaton and doesn't in- 
troduce new conserved quantities. So the parameter Q is just the temperature of the 
Schwarzschild-AdS black hole. 

2. For the U (1) ^-charged fluids discussed by [71, 72, 73] one has \& to be the bulk Maxwell 
field which gives rise to an extra parameter, q, corresponding to the black hole charge. 

3. The recent analysis of [74] incorporates three U(l) R-charges in AdSs. The matter con- 
tent of the gravitation theory comprises of three Maxwell fields and three bulk scalars. 
The black holes are described by four parameters, three charges and a temperature, 
collected here into Q in addition to the velocity field. 44 

4. One can also consider more exotic situations involving phase transitions. For instance 
in the Abelian-Higgs model coupled to gravity in AdS4 one has novel scalar hair black 
hole solutions [103]. Here *f> comprises of bulk Maxwell field and a charged scalar 
field. The set of parameters Q depends on which sector of the theory one considers 
since the bulk description involves non-trivial phase structure. In the U(l) symmetric 
phase Q just comprises of the temperature and chemical potential, whilst in the phase 
with broken U(l), one also has to incorporate the Goldstone mode. The latter phase 
gives rise to superfluid dynamics with a Landau two-fluid description. This has been 
discussed at linearized level in [104, 105] (see also [106]) and can indeed be incorporated 
into the fluid-gravity correspondence [107] (albeit with some effort since the background 
solutions are only known numerically). 

Having identified the set of parameters Q for a given example, we are in a position to 
promote them to fields depending on the coordinates x and proceed to solve the equations of 
motion (6.1) order by order in a derivative expansion. As described extensively above, this 
involves identifying two pieces of information. The first is to determine the linear operator 
H that acts on the correction terms to the metric. This operator by virtue of the Killing 
symmetries in the background spacetime (which we assume) and matter fields (6.2) is just a 
second order ordinary differential operator involving r. The second part involves identifying 
the source s n which involves a bit of work at each order. Schematically, if we promote 
Q ~^ Qi x ) an d Ufj, — > u^x) and work by expanding the solution (6.2) in the variations about 
the background value. Specifically, consider 

u, = uf + 5u,(x) = uf + e x»d u u»(x) + O (e 2 ) , 
Q = Q (0) + SQ = d v Q(x) + e a" <9 M Q (0) + O {e 2 ) . 

For simplicity of discussion consider backgrounds which are stationary, so that by passing to 
the local inertial frame we can always choose uf^ = —5^ and uf'' = x M d^fli as before using 

44 Unlike the uncharged case discussed in §5, for the fluid solutions dual to gravity which have been 
constructed so far the issue of regularity is not clear; all the analyses so far do not demand regularity of the 
Maxwell potential on the future horizon. 
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(4.7). Plugging these expansions into the metric we obtain the leading terms that contribute 
to the sources. Using 



S(r,Q(x)) =S(r,QW) + 



(oh , dS(r,Q) 



SQ (6.4) 

Q=Q(°) 



dQ 

we obtain the leading order expansion of the background metric and matter fields 

ds 2 = 2S(r,Q)dvdr-r 2 V(r,Q)dv 2 + r 2 5 ij dx i dx i 

+2 — 5Q dv dr-2S Spi dx l dr - r 2 5Q dv 2 -2r 2 <S& (1 - V)dx* dv 

* = *(r,QW) + ^6u<* + ^6Q (6.5) 

which can then be used to solve for the higher order corrections to the metric and matter 
fields This has been carried out for a wide class of models and we refer the reader to 

the original literature described at the end of §1 for details of the construction in the specific 
cases. 



Universality of transport coefficients: One of the interesting results arising from the 
explorations so far concerns the universality of transport coefficients in gravity duals. As 
is well known for two derivative gravity theories dual to boundary fluid dynamics, one has 
the famous ratio [19] of shear viscosity to entropy density saturating the conjectured bound 
v/ s ^ 1/47T. This can in fact be derived directly from the abstract discussion above, for 
the backgrounds (6.2) which respect spatial translational symmetry in global equilibrium for 
asymptotically AdS^+i spacetimes [107, 108]. 

It is interesting to ask whether the higher order transport coefficients exhibit similar 
properties. Based on the general analysis of gravity coupled to scalars and Maxwell fields 
and assuming that the theory doesn't exhibit any phase transitions (such as spontaneously 
Higgsing of the bulk gauge field), it was shown in [108] that 

A\ 1 + \ 2 = 2i 1 t^ (6.6) 

which is of course satisfied by (3.25). 



6.2 Non-relativistic fluids from gravity 

We have concentrated thus far on the dynamics of relativistic conformal fluids and con- 
structed holographic duals for them. We now briefly comment on deriving the dynamics of 
non-relativistic fluids from gravitational description. There are essentially two approaches 
taken in the literature so far: in [81] gravity duals were obtained by taking a the non- 
relativistic limit of the relativistic system, while in [80] conformal non-relativistic fluids were 
discussed in the context of Galilean holography. 
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Non-relativistic ideal hydrodynamics is described by the continuity equation, 



OtPnr + di (p nr u<) = 0, (6.7) 

together with the equation of momentum conservation, (here % = 1, . . . , d) 

d t (Pnr v l ) + <9,IF = 0, IF = p nr v* v j + 8* P nr , (6.8) 
and the equation of energy conservation, 

d t (e m + 1 p nr v 2 j + di f e = 0, f e = i (e nr + P nr ) ^ 2 . (6.9) 

where v 2 = v z v % . Here p nr is the particle number density and P nr , e nr is the pressure and 
energy density of the non-relativistic system under consideration and we use v % to denote the 
non-relativistic velocity. These equations can be derived from the relativistic stress tensor 
(2.2) by writing the conservation equation (2.1) in light-cone coordinates and demanding that 
the fluid variables be independent of one of the light-cone directions. Using the light-cone 
version of (2.1) 

d + T ++ + diT +i = , d+T +t + dfl** = , d + T + ~ + d{T^ = 0, (6.10) 

we can map the relativistic system in d dimensions into the non-relativistic equations in d — 2 
spatial directions with the following identifications: T ++ is identified with the mass density, 
T +l with the mass flux (which is equal to the momentum density), T u with the stress tensor, 
T + ~ with the energy density, and T~ l with the energy flux, 

T ++ = p nr , T +l = p nr v\ T^ = IF, 

T+- = e nr + ^ PnT v 2 , T^ = f e . (6.11) 

It is now easy to convince oneself based on (6.11) that the precise mapping between 
relativistic and non-relativistic hydrodynamic variables is 



Pnr i + i 

i / -x —p- , u =U^V , 

P = P nr , p = 2e nr + P nr . (6.12) 

The component of the relativistic velocity u~ can be determined using the normalization 
condition u^u^ = — 1 to be 

u - = ~ (-L +u + v 2 j . (6.13) 

This maps makes it clear that the transport coefficients of the non-relativistic theory 
are inherited from the parent relativistic hydrodynamics. In this description it is clear that 
non-relativistic fluids with holographic duals will saturate the conjectured viscosity bound 
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rj/s — 1/47T, which was verified explicitly in [109, 110]. Furthermore, it is also possible to 
use this light-cone reduction to infer the heat conductivity of the non-relativistic fluid: 

n nr = 2r] nr 6ni + ^ nr . (6.14) 

Pnr i 

which can be rephased as the statement that the Prandtl number of the fluid is unity. We 
recall that the Prandtl number is defined as the ratio of the kinematic viscosity u nT and the 
thermal diffusivity Xm, 

Pr = ^, (6.15) 

Xnr 

where 

Vnr = , Xnr = , (6.16) 

Pnr Pnr Cp 

where c p is the specific heat at constant pressure. 

The general idea of using slow variation in certain directions to construct inhomogeneous 
dynamical black hole spacetimes is of course not restricted to asymptotically AdS spacetimes 
and one can exploit this scheme to construct new solutions with different asymptotics. For 
general asymptotics one does not recover any interesting boundary dynamics. However, in 
recent years we have seen interesting generalizations of the AdS/CFT correspondence, such 
as the holographic duals for field theories with non-relativistic Schrodinger symmetry as 
originally discussed in [111, 112] and later embedded into string theory in [109, 110, 113]. In 
these examples one has spacetimes with non-trivial asymptotic structure and the dual field 
theory is a non-local deformation of a relativistic superconformal field theory such as M = 4 
SYM. 

In [80] the holographic duals for fluids with non-relativistic conformal symmetry i.e., 
Schrodinger invariance was constructed. In this case one can achieve the holographic dual in 
two equivalent ways: either by implementing the general procedure outlined in the lectures 
(taking into account the reduced symmetries of the problem) or equivalently by exploiting a 
specific duality transformation in string theory, the so called Null Melvin Twist [114] or the 
TsT transformation [113]. To be specific, [80] constructed fluid dynamical solutions for the 
five-dimensional effective action involving a scalar field <fi and a massive vector Am 



S = / d 5 xV^G (r - \(d M <t>)(d M <P) - -e~^F MN F MN -AA M A M - V(<j>) 

16 ttG m J V "J 4 / 

(6.17) 



where 



V(<f>) =4e 2 ^ 3 (e 2 *-4) , (6.18) 
which asymptote to the Schr 5 spacetime, 

ds 2 = r 2 (-2 dx + dx~ - (3 2 r 2 {dx + f + rfx 2 ) + (6.19) 

Consider the general form of the fluid dynamical metrics in AdSs given by (5.1) with 
d = 4. We introduce light cone coordinates x on the boundary directions. Restricting to 



48 



configurations which have cL as an isometry, after a TsT transformation on the metric (5.1) 
we get a new metric of the form 



dsl = e -3* [ -2u u Sdx fl dr + 



•>E ~ ° I A* 



Xab - — 

1 + X— 



dX A dX B 



1 + X- 



(6.20) 



with 



XA = XA-~u.S6 r A . (6.21) 

The TsT transform converts the asymptotically AdSs spacetime (5.1) to an asymptotically 
Schr 5 spacetime, which depends on the paramters b, j3, Ui which are arbitrary functions of 
(x + ,x). The relativistic velocity field w M in fact descends naturally into a non-relativistic 
velocity field Vi via the map 



u + = /3, u { = (3 



(6.22) 



which can be inferred from the light-cone reduction of the stress tensor [80]. In the field 
theory description b is related to the temperture, while (3 sets the chemical potential for the 
non-relativistic particle number. 

Conformal non-relativistic fluids with Schrodinger symmetry are in general compressible; 
this follows from the fact that the energy density is related to the pressure 2e nr = dP nr 
through the equation of state (which in turn follows from scale invariance). To make contact 
with the usual studies of incompressible Navier-Stokes equations we need to ensure that we 
can decouple the fluctuations in the density. This can be achieved by looking at low frequency 
modes which do not excite the propagating sound mode in a hydrodynamic system, i.e., by 
focussing on the shear mode. In fact, this limit was discussed recently in the context of the 
fluid-gravity correspondence in [81] (see also [115, 116] for closely related results), where the 
authors showed that starting from a parent relativistic conformal fluid dynamical system 
one can recover incompressible Navier-Stokes equations in a suitable scaling limit. This 
scaling limit reveals an interesting structure in the fluid equations - they are scale invariant 
under a new scaling symmetry, which is the Galilean conformal algebra discussed recently in 
[117]. This symmetry is different from the Schrodinger symmetry enjoyed by the conformal 
non-relativistic fluids discussed above (see also [118], [119]). 

One can view these two constructions as follows: given a relativistic theory, in particular 
relativistic fluid dynamics, one can attain a non-relativistic limit either by (i) taking the 
speed of light to infinity or (ii) by a light-cone reduction. The former procedure is related 
to contracting the Poincare algebra to the Galilean algebra and when implemented on the 
relativistic Navier-Stokes equations leads to the incompressible non-relativistic Navier-Stokes 
fluid of [81]. The latter procedure of light-cone reduction converts d dimensional relativistic 
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fluid dynamics into a d — 2 spatial dimensional non-relativistic fluid dynamics. Since one 
requires a null vector to reduce on the light-cone, we end up losing an effective dimension in 
our field theoretic description. 

7 Discussion 

In the course of these lectures we have discussed the essential features which relate the 
physics of inhomogeneous, dynamical, black hole solutions of general relativity in asymptot- 
ically AdS spacetimes, to fluid dynamics of strongly coupled superconformal field theories 
living on the boundary of these AdS geometries. In particular, given any solution to the 
relativistic Navier-Stokes equations and their non-linear generalizations (with the strong cou- 
pling values of the transport coefficients), one can construct an asymptotically AdS black 
hole solution. Alternately, the general construction described here can be viewed as a deriva- 
tion of the hydrodynamic stress tensor of the dual superconformal field theories. While we 
have discussed the construction to second order in derivatives, it is clear that the general con- 
struction can in principle be extended to arbitrary orders in the gradient expansion (albeit 
with increasing computational complexity in evaluating the source terms) . 

We have also described how one can use the geometric description to understand the field 
theory entropy and constructed a specific gravitational entropy current. An important issue 
which we have briefly discussed concerns the geometric description of entropy and subtleties 
associated with identifying the entropy with the area of the event horizon. While it seems 
reasonable to associate the area of the event horizon to the field theory entropy in situations 
where one has slow variations, it seems clear that this cannot be true in general as exemplified 
by the conformal soliton flow described in [63]. 

Furthermore, we have reviewed various generalizations of the correspondence over the 
past year or so which have led to interesting new insights into forced fluid dynamics and 
charge transport, etc.. These provide an interesting arena for further exploration - in princi- 
ple it should be possible to derive the complete stress tensor and charge currents for M = 4 
SYM up to two derivatives incorporating the three U(l) R-charge and angular momentum 
chemical potential. In addition, it is clear that there are many interesting directions that 
can be tackled within this framework, most notably the issues which were raised in §1 as 
part of the motivation for the correspondence. 

One issue we have not touched upon in these lectures is the relation of the fluid-gravity 
correspondence to the membrane paradigm [9]. Both purport to identify the dynamics of 
black holes with hydrodynamics, and indeed the correspondence offers a important new per- 
spective on this issue [120]. 45 Since in the fluid-gravity correspondence the entire spacetime 
dynamics is mapped unambiguously into the boundary fluid dynamics, it is natural to argue 
that the correspondence in fact implements the ideas inherent in the membrane paradigm, 

45 An earlier discussion of the membrane paradigm in the context of hydrodynamics in AdS/CFT can be 
found in [33]. Also sec [121] for a recent discussion of the membrane paradigm in the holographic context. 
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albeit with a new wrinkle: the membrane is not located at the stretched horizon, but rather 
at the boundary of the spacetime. Despite its location, this 'membrane at the end of the 
universe' faithfully captures the entire bulk spacetime dynamics and implements the mem- 
brane paradigm holographically, as is clear from the fact that we not only recover the stress 
tensor of the fluid, but also a particular gravitational entropy current. Note that in the 
conventional membrane paradigm, the dynamics of the stretched horizon only captures the 
physics of the region behind the black hole horizon and not of the entire spacetime. One 
can of course use the AdS/CFT correspondence to interpolate between these two extreme 
descriptions [122]: if one so wishes, it is possible to define a fiducial membrane at some other 
radial location (and in particular on the stretched horizon) and identify the fluid dynamics 
on this surface. The dynamics on this surface is related to the boundary fluid dynamics by 
an appropriate flow equation which can be derived from the bulk gravitation equations using 
techniques analogous to the holographic renormalization group. 

Finally, we should mention recent investigations of gravitational duals of field theories 
perturbed away from equilibrium [101, 123]. As we have emphasized the hydrodynamic 
behaviour comes into play only when the field theory achieves local thermal equilibrium. 
In general, large perturbations will evolve outside the fluid regime, till such a time the 
system thermalizes and achieves local equilibrium. Correspondingly, in the dual gravity 
description, black hole formation is outside the long-wavelength hydrodynamic regime. The 
remarkable aspect of the recent discussion of [123] is that the system tends to "thermalize" 
rapidly - in effect the hydrodynamic description takes over almost instantaneously after the 
perturbation is switched off. 46 It would be interesting to understand this passage to local 
thermal equilibrium for a generic perturbation in a strongly coupled field theory in greater 
detail. 
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A Weyl covariant curvature tensors 



One can define a curvature associated with the Weyl covariant derivative by the usual pro- 
cedure of evaluating the commutator between two covariant derivatives. Defining a field 
strength for the Weyl connection 

J> = V M A - VuA» (A.l) 
we find a Riemann curvature tensor IZ^xa' 

TZ^uXa = R^v\a + 4 ^g u][x S^ \V a Ap + A a Ap - ^g a p\ - J> g\ a (A. 2) 

These two tensors are in an appropriate sense Weyl invariant; it is not hard to check from 
the definitions that J 7 ^ = T^y and TZ^ u \ a = TZ^ u x a - It should be borne in mind that the 
Weyl covariant Riemann tensor defined above (A. 2) has different symmetry properties from 
the conventional Riemann tensor. Most importantly, it is not anti-symmetric under the 
exchange of the last two indices. For example, 

TZfivXa T^fivcrX 2.7"" ul/ Q\a 

TZ-fiuXa — TZ-Xafiu = ^9u][X^] ~ F uu 9Xu + F\o 9fiu (A. 3) 

n, auP v a - n ua ^ v a v 13 = -t^ v a v a 

With the basic definitions (A.l) and (A. 2) in place we can proceed to define analogous 
expressions for various other tensors often encountered in differential geometry, such as the 
Ricci tensor, Ricci scalar: 

TZfiu — TZfiav Rfiv (d 2) (V fj,A v + A^A V A g^v} g^v^ xA ■F^v T^\w , . , 
K = TZ a a = R - 2(d - 1)V X A X + (d-2)(d- 1)A 2 = e~ 2 ^TZ . 

In addition it is also worth noting that the conformal tensors of the underlying spacetime 
manifold appear as a subset of conformal observables in hydrodynamics. These conformal 
tensors are the Weyl-covariant tensors that are independent of the background fluid velocity, 
for we have already accounted for the terms involving velocity derivatives explicity above. A 
well known example of this class of operators is the Weyl curvature C^xa which is the trace 
free part of the Riemann tensor. In d > 3 it is given as 

CfMuXa = Tl^uXa + ^ ^9v][x8 a ^S a p = C^Xa ~ F^u9Xa = CfivXcr (A-5) 

where we introduced the Schouten tensor S^ v is defined as 47 
<V = ~T~~~n Up, - = S, v - (v,A v + A,A V - ^) -P^=S, V (A.6) 



d-2 v " 2(d-i)y ^ yr**,-^^ d _ 2 
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From equation (A. 5), it is clear that C^\ a = C^xa + Tp, v gxcj is clearly a conformal 
tensor. Such an analysis can in principle be repeated for the other known conformal tensors 
in arbitrary dimensions. 

The Weyl Tensor C^xa has the same symmetry properties as that of Riemann Tensor 

C^i/Xu Cjy^Xu C^i/ a x Cxafiu . . 

(A. 8) 

and C m x a = . 

It follows that C^ av p u a viP is a symmetric traceless and transverse tensor, which is why it 
shows up as a second order contribution in the conformal hydrodynamical stress tensor. 
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